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ABSTRACT 


Analyses  and  measurements  of  attenuation  of  thermal  radia¬ 
tion  leaving  a  heated  plate  caused  by  a  layer  of  absorbing  gas 
are  discussed.  Both  a  near  stagnant  fluid  and  flow  over  a  flat 
plate  were  investigated  by  analysis  and  experiment.  An  appara¬ 
tus  in  which  various  porous  plates  were  heated  electrically  was 
designed  and  constructed  so  that  it  could  be  adopted  to  either 
the  stagnant  or  flow  systems.  Radiant  energy  fluxes  in  speci¬ 
fied  wavelength  bands  leaving  the  plate  were  measured  with  an 
infrared  radiometer  developed  and  built  for  this  purpose. 

Specifically,  radiometric  measurements  for  the  4.3y  band 
of  carbon  dioxide  were  made  for  various  plate  surface  tempera¬ 
tures,  injection  rates,  thicknesses  of  the  absorbing  gas  layer 
and  viewing  angles.  In  addition  to  the  above  parameters  the 
Reynolds  number  was  varied  and  the  effect  of  the  starting  length 
investigated  for  flow  over  a  flat  plate.  It  was  found  that  the 
gas  is  capable  of  attenuating  a  significant  fraction  of  the 
radiation  leaving  the  hot  plate.  Increases  in  the  plate  tem¬ 
perature  and  gas  layer  thickness  increased  the  attenuation. 

For  a  given  mass  injection  rate  of  an  absorbing  gas  the  attenu¬ 
ation  increased  with  an  increase  with  the  Reynolds  number.  The 
analyses  and  experiments  show  that  in  order  to  obtain  signifi¬ 
cant  attenuation  over  large  spectral  regions  a  mixture  of 
selected  gases  or  gas  particulate  mixtures  would  have  to  be 
injected. 
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1.  INTRODUCTION 


The  emission  of  thermal  radiation  from  high  temperature 
bodies  has  become  of  concern  to  design  engineers.  The  struc¬ 
tural  components  of  a  jet  engine,  for  example,  and  the  exhaust 
gases  emit  not  only  radiation  in  the  visible  but  also  in  the 
infrared  range  of  the  spectrum.  Furthermore,  the  soot  and 
carbon  particles  suspended  in  the  exhaust  stream  contribute 
to  the  radiation.  The  shielding  of  components  and/or  the 
environment  from  intense  radiation  can  be  an  important  objective 
in  the  design  of  aircraft  and  missiles,  airfoil  surfaces,  etc. 

In  most  applications,  shielding  is  used  to  reduce  the  heating 
or  cooling  of  a  system  component;  but  radiant  energy  received 
by  a  sensor  from  a  source  may  also  be  reduced  by  similar  means 
to  decrease  the  probability  of  detection. 

The  infrared  radiation  leaving  a  source  may  be  decreased, 
for  instance,  by  geometrical  design  and/or  reduction  of  emit- 
tance.  The  radiation  can  also  be  reduced  by  transpiration  of 
radiation— absorbing  gas  into  the  gas  stream  through  a  porous 
material.  The  injection  of  an  absorbing  fluid  into  the  gas 
stream  may  shield  the  hot  surfaces  from  an  observer  and  will 
also  cool  the  surfaces  and  thus  reduce  the  radiant  energy 
emitted.  The  extent  of  the  second  benefit  would  be  increased 
by  injecting  a  cold  liquid  which  would  vaporize  as  it  transpires 
through  the  porous  structure.  A  liquid  additive  might  be  pre¬ 
ferred  over  a  gaseous  one  since  the  storage  volume  for  a  given 
mass  would  be  an  order  of  magnitude  smaller.  The  infrared 
suppression  scheme  would  be  used  only  when  need  arises,  and 
the  consumption  of  the  radiation  absorbing  fluid  could  be  kept 
at  a  minimum. 

The  present  study  has  been  conceived  as  an  analytical  and 
experimental  investigation  of  the  feasibility  of  attenuating 
the  infrared  radiation  emitted  from  hot  surfaces  by  injection 
of  an  absorbing  gas  between  the  source  and  the  detector.  The 
objective  of  the  research  program  was  to  determine  how  effec¬ 
tively  the  transpiration  of  an  absorbing  gas  attenuates  the 
radiation.  The  purpose  of  the  investigation  conducted  was  two¬ 
fold  and  can  best  be  described  by  briefly  summarizing  the 
essential  elements  of  the  research  program: 

1.  Prediction  of  attenuation  of  radiation  emitted  from 
a  hot  surface  through  which  an  absorbing  emitting 
gas  is  injected. 
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2.  Performance  of  experiments  to  verify  the  analytically 
predicted  attenuation. 

The  program  as  stated  above  is  too  ambitious  and  too 
inclusive  to  be  treated  thoroughly  within  the  limits  of  the 
time  and  funding  of  the  program.  The  specific  tasks  under¬ 
taken  were  the  following: 

1.  Prediction  and  measurement  of  attenuation  by  a  near- 
stagnant  layer  of  absorbing  gas. 

2.  Prediction  and  measurement  of  attenuation  of  infrared 
radiation  emitted  by  a  hot,  flat  plate  through  or 
ahead  of  which  an  absorbing  gas  is  injected  into  the 
stream. 

In  Chapter  2  the  problem  is  discussed  in  more  detail  and 
the  literature  is  reviewed.  Chapter  3  describes  the  analysis 
and  apparatus  used  for  the  investigation  of  attenuation  of 
radiation  by  a  near-stagnant  layer  of  radiating  gas  and  the 
comparison  of  the  results  is  given.  In  Chapter  4,  an  analysis 
is  presented  and  the  experimental  facility  is  described  for 
the  study  of  attenuation  of  radiation  in  boundary  layer  flow 
over  a  flat  plate.  Detailed  comparison  is  made  between  pre¬ 
dicted  attenuation  and  experimentally  obtained  results. 

Chapter  5  contains  conclusions  and  recommendations. 
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2. 


DISCUSSION  OF  THE  PROBLEM 


2 . 1  Attenuation  of  Infrared  Radiation  by  Gases 

To  be  specific,  consider  a  plane  layer  of  absorbing- 
emitting  gas  of  thickness  6  which  is  in  contact  with  a  hot 
surface.  The  spectral  intensity  of  radiation  emerging  from 
the  top  of  the  layer,  Ix(6,8),  neglecting  scattering  and^ 
assuming  local  thermodynamic  equilibrium  is  given  by  [1] 

6  6 
-/„  K.dy/cose  r6  -I  ic.dt/cose 

Ix(6,0)  =  Iu(0,6)e  ♦  !bx(y)e  Y  Kxdy/cos0 

*  0 

(2.1) 


In  writing  this  equation  we  have  assumed  that  the  surface  and 
the  layer  are  axisymmetric ,  i.e.,  the  radiation  surface  charac¬ 
teristics  and  the  spectral  intensity  are  independent  of  the 
azimuthal  angle.  The  physical  meaning  of  Eq.  (2.1)  is  clear: 

It  expresses  the  fact  that  the  intensity  emerging  from  the 
layer  in  a  given  direction,  0,  results  from  a  contribution  due 
to  radiation  leaving  the  boundary  plane,  y  =  0,  attenuated  by 

the  factor  exp ( -/O^xxdy/cos0]  and  contributions  from  gaseous 


emission  at  all  interior  points,  y,  reduced  by  the  factor 
exp[-/  tc,dt/cos0]  to  allow  for  absorption  by  the  intervening 

y  A 

gas.  The  spectral  intensity  of  radiation  leaving  the  surface 
IU(O,0)  is  the  sum  of  the  emitted  and  reflected  radiation  and 
can  be  expressed  as  (1] 


iu«M).cx(e)ibx(y2j 
f6 

■  I  'bx^ 


ir/z 


'ex'9’' 


/,  *xdy/cos0' 


(2.2) 


fQ  <xdt/cos0' 


<xdy/cos0' 


Px  (O^OjsinO ' cos6 'd0 ' 


where  ex(8)  is  the  directional  emittance,  px(8‘-*8)  is  the  bi¬ 
directional  reflectance,  and  Ie^  is  the  intensity  of  external 
incident  radiation.  The  terms  in  Eq.  (2.2)  represent  the 
surface  emission,  the  reflection  of  the  energy  incident  on  the 


*  Numbers  in  brackets  denote  references. 
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layer  after  being  attenuated  by  the  layer  and  the  reflection 
of  attenuated  gaseous  emission,  respectively.  For  a  trans¬ 
parent  layer  of  gas  <A  =  0  and  Eq.  (2.1)  reduces  to 

Ia(6,0)  =  ln(0,e)  (2.3) 

The  effectiveness  of  a  layer  of  a  participating  gas 
placed  over  a  surface  to  absorb  radiation  can  be  expressed 
by  the  normalized  spectral  intensity,  which  is  defined  as  the 
ratio  of  intensity  at  the  edge  of  the  layer  to  the  intensity 
leaving  the  surface  if  the  gas  layer  was  transparent,  i.e., 


lj(fi,e)  =  Ix(6,0)/In(O,0) 


Transparent 


(2.4) 


and  serves  as  a  measure  of  the  spectral  attenuation.  Gases 
absorb  and  emit  radiation  in  bands,  and  hence  the  normalized 
intensity  for  band  i  is  given  by 

I*  (5,0)  »  f  I* (5 , 6) dX  (2.5) 

l  JAX^ 


The  linear  (volumetric)  spectral  absorption  coefficient, 
■c.,  does  not  only  depend  on  the  nature  of  the  gas  but  on  con¬ 
centration,  temperature  and  total  pressure  of  the  gas.  The 
spectral  absorption  coefficient  can  also  be  expressed  as 

K\  a  *  picA  (2-6) 

where  *A  and  tc'A'  are  referred  to  as  the  spectral  mass  absorption 
coefficients.  The  quantity  pxd/Ng,  where  Ng  is  the  number  den¬ 
sity  (number  per  unit  volume)  “of  gas  molecules,  has  the  dimen¬ 
sion  of  an  area  and  can  be  interpreted  as  a  cross-section  for 
absorption  of  radiation  of  wavelength  X.  Since  the  Planck's 
function  I^x  is  a  function  of  temperature  it  is  clear  from 
Eqs.  (2.4),  (2.5)  and  (2.6)  that  the  prediction  of  the  normal¬ 
ized  (spectral  or  band)  intensity  emerging  from  the  gas  layer 
is  straight  forward  once  the  heat  and  mass  transfer  problem 
has  been  solved. 

2. 2  Absorption  Characteristics  of  Gases 

The  description  of  the  radiation  properties  of  gases  from 
a  quantum  mechanical  viewpoint  is  discussed  in  Reference  (2]. 
The  problem  of  predicting  the  properties  of  gases  is  very  com¬ 
plex  and  has  been  the  subject  of  much  recent  work  {2,  3].  The 
experimental  evaluation  of  the  spectral  absorption  coefficient 
is  extremely  difficult  and  time  consuming  because  of  a  large 
number  of  independent  parameters.  In  theory,  it  is  possible 
to  predict  it  from  quantum  mechanics;  however,  for  gases  the 
complexity  of  the  calculations  makes  this  approach  in  general 
impractical.  Hence  the  absorption  properties  are  conveniently 
taken  from  experimental  data  and/or  from  analytical  models  (3). 
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An  extensive  bibliography  on  infrared  radiative  properties  of 
gases  has  been  compiled  recently  [4],  and  there  is  no  need  to 
repeat  that  survey. 

A  literature  search  was  undertaken  to  compile  a  list  of 
gases  which  have  desirable  absorption  characteristics  in  the 
infrared  region  of  the  spectrum  and  could  be  employed  as  in¬ 
jected  gases.  The  common  gases  such  as  oxygen,  nitrogen,  and 
hydrogen  have  symmetrical  molecules  and  for  all  practical 
purposes  are  transparent  to  thermal  radiation.  On  the  other 
hand,  gases  and  vapors  such  as  carbon  dioxide,  sulfur  dioxide, 
ammonia,  water  vapor,  carbon  monoxide  and  various  hydro¬ 
carbons  absorb  and  emit  thermal  radiation  in  bands.  For 
example,  the  important  radiation  bands  for  C02,  lie  in  wave¬ 
length  region,  2.36  to  3.02y,  4.01  to  4.80y,  and  12.5  to 
16. 5y  [5j.  A  list  of  28  absorbing  gases  is  presented  in  Table 
1.  Figures  1  through  7  illustrate  the  infrared  transmission 
spectra  of  these  gases.  All  of  the  spectral  data  were  obtained 
with  a  Perkin-Elmer  Model  21  spectrometer  having  the  following 
features;  sodium  chloride  optics,  10-cm  gas  layer  cell  at  room 
temperature  (23°C)  with  sodium  chloride  windows,  and  a  wave¬ 
length  range  from  2  to  15  y  [6j. 

Table  2  contains  brief  qualitative  remarks  as  to  the 
toxicity,  fire,  and  explosive  hazards  of  the  gases.  The  infor¬ 
mation  contained  in  Table  2  was  taken  from  Reference  [7],  It 
should  be  noted  that,  in  general,  those  substances  classified 
as  having  "slight  toxicity"  produce  changes  in  the  human  body 
which  are  readily  reversible  and  which  will  disappear  following 
termination  of  exposure,  either  with  or  without  medical  treat¬ 
ment.  Those  substances  classified  as  having  "moderate  toxicity" 
may  produce  irreversible  as  well  as  reversible  changes  in  the 
human  body.  These  changes  are  not  of  such  severity  as  to 
threaten  life  or  produce  serious  permanent  physical  impairment. 
"Severe  toxicity"  includes  those  substances  which  on  a  single 
exposure  lasting  seconds  or  minutes  can  cause  injury  sufficient 
to  threaten  life  or  to  cause  permanent  impairment,  disfigurement, 
or  irreversible  changes. 

Since  the  absorption  properties  of  carbon  dioxide  are  well 
known  and  the  gas  is  relatively  non-hazardous ,  it  was  chosen  as 
working  fluid  for  this  investigation.  The  fundamental  4.3y 
band  was  chosen  for  the  study.  The  choice  of  the  gas  and  the 
band  was  somewhat  arbitrary. 

2 . 3  Heat  Transfer  in  Radiating  Gases 

The  presence  of  a  finite  interfacial  velocity  normal  to 
a  surface  where  simultaneous  heat,  mass  and  momentum  transfer 
are  taking  place  alters  appreciably  the  temperature,  concentra¬ 
tion  and  velocity  profiles.  It  has  been  found  that  one  of  the 
most  effective  ways  of  cooling  a  high  performance  surface  is 
that  of  mass  transfer  cooling  which  involves  the  injection  of 
a  foreign  gas  (or  liquid)  into  the  boundary  layer.  A  detailed 
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Table  1.  Spectral  Absorption  Characteristics 
of  Various  Gases 


Spectrum 

No. 

Compound 

Chemical  Formula 

3 

Acetylene 

HC  :  CH 

S 

Allene 

CH2  :  C  :  CH2 

6 

Ammonia 

NHj 

9 

1.3  Butadiene 

CH2  :  CH  •  CH  :  CH2 

10 

n- Butane 

CH,CH2CH2 CHj 

11 

Carbon  Dioxide 

C02 

12 

Carbon  Monoxide 

CO 

13 

Cyanogen 

(CN)2 

14 

Cyclopropane 

(ch2), 

15 

Diazomethane 

ch2n2 

18 

Dimethylamine 

(CHs)2NH 

19 

Dimethyl  Ether 

CHjO  CH, 

22 

Ethane 

CH,CH, 

28 

Ethylene 

CH2  :  CH2 

38 

Isobutane 

(CH,),  CH 

39 

Methane 

CHS 

42 

Methyl  Chloride 

CH,C1 

43 

Methyl  Mercaptan 

CH,SH 

45 

Monoethylamine 

c2h5nh2 

46 

Monomethy lamine 

ch,nh2 

48 

Nitric  Oxide 

NO 

56 

Nitrous  Oxide 

N,0 

58 

Propane 

ch,ch2ch, 

59 

Propylene 

CH,  •  CH  :  CH, 

60 

Propyne 

CH,  •  C  :  CH 

61 

Sulfur  Dioxide 

SO, 

64 

Trine thylamine 

(CH,),N 

66 

Vinyl  Chloride 

CH,  :  CH  .  Cl 

Spectral  Absorptance  of  Acetylene,  AUene 
Aaaonia,  and  I  .  J  butadiene 


Um  WiUf 

mm. 

Rr*4i«lk« 

I'iguro 


£  "Har  b  "HA  i  at  i  "HA  &  & 

3.  Spectral  Absorptancc  of  Cyclopropane,  Uiazo- 
methane,  Diucthyl  Amine,  and  Dimethyl  Htner 


Spectral  Ahsorpt  once  «»l  Mothvl  Chloride.  Methyl 
Mercaptan,  Monoethyl  Amine,  and  Monometnyl  Amine 
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Figure  7.  Spectral  Absorptance  of  I’ropyne,  Sullur 

Dioxide,  Tr  imetliv  lamme  ,  and  Vinyl  Chloride 


Table  2.  Hazardous  Properties  of  Various  Gases 


Gas 

Toxicity 

Fire 

Hazard 

Explosive 

Hazard 

Acetylene 

Slight-Moderate 

Dangerous 

Moderate 

Allene 

Unknown 

Dangerous 

Moderate 

Ammonia 

Severe 

Moderate 

Moderate 

1.3  Butadiene 

Moderate 

Dangerous 

Moderate 

n- Butane 

Slight-Moderate 

Dangerous 

Moderate 

Carbon  Dioxide 

Slight 

Non 

Flammable 

- 

Carbon  Monoxide 

Severe 

Dangerous 

Severe 

Cyanogen 

*  Severe 

- 

- 

Cyclopropane 

Moderate 

Dangerous 

Moderate 

Diazomethane 

Moderate 

- 

Severe 

Dimethylamine 

Moderate 

Dangerous 

Moderate 

Dimethyl  Ether 

Moderate 

Highly 

Dangerous 

Moderate 

Ethane 

Slight-Moderate 

Dangerous 

Moderate 

Ethylene 

Slight-Moderate 

Dangerous 

Moderate 

Isobutane 

Slight 

Dangerous 

Severe 

Methane 

Slight 

Dangerous 

Dangerous 

Methyl  Chloride 

Severe 

Dangerous 

Moderate 

Methyl  Mercaptan 

Moderate 

Dangerous 

- 

Monoethylamine 

Unknown 

Unknown 

- 

Monome  thy 1  amine 

Severe 

Dangerous 

Moderate 

Nitric  Oxide 

Severe 

Unknown 

- 

Nitrous  Oxide 

Slight-Moderate 

Moderate 

Moderate 

Propane 

Slight 

Highly 

Dangerous 

Severe 

Propylene 

Slight-Moderate 

Dangerous 

Moderate 

Propyne 

Slight-Moderate 

Dangerous 

Moderate 

Sulfur  Dioxide 

Severe 

Non 

Flammable 

- 

Trimethylamine 

Moderate 

Moderate 

Moderate 

Vinyl  Chloride 

Sevore 

Dangerous 

Severe 
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description  concerning  specific  examples  of  some  mass  transfer 
cooling  processes  can  be  found  in  Reference  [8].  A  method  of 
particular  interest  in  this  investigation  is  that  of  transpira¬ 
tion  cooling  associated  with  the  injection  of  coolant  gas 
through  a  porous  structure  into  the  surrounding  high  tempera¬ 
ture  boundary  layer.  A  very  extensive  (over  100  references) 
bibliography  of  studies  relating  to  transpiration  cooling  has 
been  compiled  recently  by  Nealy  [9]  and  will  not  be  repeated 
here;  instead  the  emphasis  will  be  placed  on  the  interaction 
of  convection  and  radiation  when  the  injected  gas  is  capable 
of  absorbing  and  emitting  radiation. 

The  problem  of  convection-radiation  interaction  has  been 
a  subject  of  considerable  study  during  the  last  decade.  A 
summary  of  these  studies  dealing  primarily  with  external  type 
flow  of  gray  gases  without  mass  injection  can  be  found  in 
publications  [1,  4,  10].  The  Couette  flow  geometry  has  re¬ 
ceived  considerable  attention  [11-13]  primarily  because  it 
represents  a  simple  model  of  boundary  layer  flow.  The  analyti¬ 
cal  studies  [14-17]  dealing  with  the  radiation  effects  in 
boundary  layer  flow  over  a  flat  plate  of  an  absorbing  gas  have 
been  extended  by  Kennedy  [18]  to  account  for  mass  addition. 

The  effect  on  the  convection  and  radiation  heat  transfer  by 
injection  of  an  absorbing  gas  in  two-dimensional  and  axisym- 
metric  stagnation  point  regions  has  been  reported  by  Macken 
and  Hartnett  [19].  Howe  [20,  21]  with  a  specific  problem  in 
mind  studies  the  effect  of  foreign  species  on  the  convection- 
radiation  for  an  axisymmetric  stagnation  region  in  the  presence 
of  a  shock  front;  his  results  are  applicable  only  to  his 
specific  case  studied. 

In  the  above  references  the  gray  gas  assumption  has  been 
employed.  Although  physically  questionable,  the  simplifications 
have  been  necessitated  by  the  lack  of  spectral  data  for  the 
absorption  coefficients  and  the  need  to  reduce  the  complexity 
of  the  conservation  equations.  The  utility  of  this  approach 
has  been  that  the  simplified  solutions  illustrated  the  effects 
of  various  parameters.  Only  very  recently  has  the  gray  assump¬ 
tion  been  relaxed  and  a  number  of  solutions  reported  [22-26], 

In  all  of  these  studies  the  fluid  was  considered  to  be  stagnant 
and  use  was  made  of  the  total  band  absorptance  correlation  for 
the  gases  [3].  The  solutions  reported  in  References  [25,  26] 
were  for  linearized  radiative  transfer.  Nongray  analyses  of 
radiant  heat  transfer  for  flow  of  high  temperature  air  over  a 
flat  plate  without  mass  injection  [27]  and  over  symmetric 
bodies  with  massive  blowing  [28,  29]  have  been  reported. 

To  the  authors'  knowledge  no  pertinent  experimental  studies 
dealing  with  combined  convection  and  radiation  heat  transfer 
with  transpiration  of  an  absorbing  gas  have  been  published. 

Some  related  experiments  dealing  with  combined  conduction  and 
radiation  through  a  stagnant  plane  layer  of  infrared  radiation 
absorbing  gas  confined  between  two  opaque  walls  have  been 
carried  out  [30,  31];  and  the  analysis  based  on  wide  band  total 
absorptance  model  have  been  substantiated. 
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3.  ATTENUATION  OF  RADIATION 
BY  NEAR-STAGNANT  ABSORBING  GAS 


3. 1  Analysis 

3.1.1  Physical  Model  and  Assumptions 

The  physical  model  together  with  its  coordinate  system  is 
shown  in  Figure  8.  An  absorbing- emitting  gas  is  injected 
steadily  and  uniformly  downward  through  a  heated  porous  plate. 
The  hot  gas  is  held  adjacent  to  the  plate  by  buoyancy  forces 
aided  by  retaining  sidewalls  which  are  assumed  not  to  partici¬ 
pate  in  the  radiant  energy  transfer.  The  plate  dimensions  are 
large  compared  to  the  layer  thickness,  so  that  the  flow  and 
heat  transfer  can  be  considered  to  be  one-dimensional.  The 
thin  two-dimensional  layer  in  which  the  gas  escapes  tne  side- 
walls  is  neglected.  This  neglect  is  justifiable  as  long  as 
the  injection  rates  are  small.  Surrounding  the  system  is  a 
medium  which  is  transparent  to  radiation  of  all  wavelengths 
investigated  here.  A  spectral  radiation  intensity  originating 
from  sources  external  to  the  system  and  incident  of  the  gas 
layer  is  denoted  by  Ie(j).  Only  binary  mixtures  composed  of  the 
injected  and  surrounding  gases  are  considered.  The  mechanisms 
of  diffusion  are  ordinary  and  thermal  diffusion;  pressure  and 
body- force  diffusion  are  neglected. 


3.1.2  Conservation  Equations 


In  the  steady  state,  conservation  of  mass  is  satisfied  by 
a  constant  mass  flux  through  the  layer  equal  to  the  injected 
mass  flux.  The  momentum  conservation  equation  for  the  physical 
model  is 

"a?  ■  '  *  gSpCT.-T).  (3.1) 

Conservation  of  species  is  expressed  by 

37  [""i  *  h]  ■  0  (3-2) 

where  the  diffusional  mass  flux  of  the  i ^  species  is  given  by 


wid-wi)  dT 
- T -  3y.  * 


(3.3) 
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Figure  8.  Physical  Model  and  Coordinate  System 


The  second  term  in  the  above  equation  represents  the  thermo¬ 
diffusion  (T-D)  effect  [32].  The  energy  equation  as  applied 
to  the  system  is  expressed  as 


f 

- 

f 

1 

”CP  *  i > 

C  *  C  ' 

Pl  P2^ 

T  ♦  q 

J 

d 

37 


where  the  heat  flux  is  given  by 

q  =  'R3y 


=  o 


v  dT  *  r  *  DTf  M‘  ■ 

k  ^  +  F  +  RT|NT^7JaJi 


(3.4) 


(3.5) 


Included  in  the  energy  equation  are  the  convection,  inter¬ 
diffusion,  conduction,  radiant  and  di ffus ion- thermo  (D-T)  heat 
transfer  mechanism  terms. 


Boundary  conditions  at  the  edge  of  the  one-dimensional 
layer,  6,  are  taken  to  be  the  same  as  those  at  the  outer  edge 
of  the  thin  two-dimensional  layer.  Due  to  the  diffusion  and 
blowing,  the  gas  temperature  immediately  adjacent  to  the  plate 
may  differ  from  the  plate  temperature  [32,  33].  Hence  the 
boundary  conditions  for  the  conservation  equations  are 


y  -  0;  T  -  T0  -  Tr 
y  -  6;  P  *  P„,  T  -  T„,  w,  -  wloe. 


(3.6) 


Assuming  that  the  porous  plate  is  impermeable  to  the  surround¬ 
ing  species  provides  the  remaining  required  boundary  condition. 
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The  species  (3.2)  and  energy  (3,4)  equations  were  inte¬ 
grated  once  to  reduce  the  order  of  the  equations.  For  the 
species  equation,  the  resulting  integration  constant  is  derived 
from  the  impermeable  plate  condition  and  found  to  be  the  net 
injected  mass  flux,  m.  In  the  energy  equation  the  integration 
constant  is  the  energy  flux,  ew,  which  is  constant  under  steady 
state  conditions.  Defining  dimensionless  variables  as  the 
ratio  of  the  actual  variable  to  a  reference  value  and  collect¬ 
ing  terms,  the  first  integrals  of  the  conservation  relations 
take  the  following  dimensionless  forms, 


Momentum: 
Species : 
Energy :  c* 


dp  + 

ar 


+  2 

P 


P+(0-0J 


W1 


pV  dwi  A  awi ri-*i)  do 
=  W  - (5 - 1 


G  Ac  + 
P 


NM+G+ 


0 


k+  dQ 

Pe  ar 


♦  BoF+ 


(3.7) 

(3.8) 

(3.9) 


with  boundary  conditions 


S  =  0;  0  -  O0  =  1 

c  •  1;  p*  -  1,  0  »  ©„,  w j  =  wloo. 
A  parameter,  N,  is  defined  as 


(3.10) 


N  =  Ro/c  =  RaT/c  T  (3.11) 

r  r 

which  in  light  of  Eq .  (3.5),  represents  the  ratio  of  energy 
transport  due  to  the  di f fusion- thermal  effect  to  the  energy 
transport  by  convection.  Since  diffusion- thermo  is  a  low- 
order  effect,  this  parameter  becomes  significant  only  when  the 
convection  and  di f fusion- thermo  energy  transport  are  of  the 
same  order  of  magnitude  as  in  free  convection. 


To  complete  the  system  of  equations,  expressions  for  the 
radiation  flux  and  an  equation  of  state  are  required.  In 
general,  the  properties  are  functions  of  the  thermodynamic 
state.  For  simplicity,  however,  the  necessary  transport 
properties  are  considered  to  be  only  temperature  dependent  and 
are  expressed  by  simple  power  law  relations 

k*  ■  9b,  D*  »  0C  (3.12) 


where  b  and  c  are  selected  for  each  gas  over  the  temperature 
range  of  interest. 

5.1.3  Radiation  Flux  and  Emergent  Intensity 

In  the  absence  of  scattering  the  spectral  radiation  flux 
is  (I] 
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0 

with  the  optical  depth  and  thickness  being  defined  by 

fy  r6 

t  =  i  tc  dy  and  =  tc  dy  (3.14) 

u  I A  to  7  otojo)7  v  7 

J0  '0 

respectively . 

The  intensity  of  radiation  leaving  the  plate  which  is  the 
sum  of  the  emitted  and  reflected  contributions  is  given  by 
Eq.  (2.2).  In  terms  of  new  variables  it  can  be  written  as 

=  StfMWy  +  2j  iP(i)(u,^)Ieti)(u,)exp(T0a)/u-)u'dii' 

■2  j  |J  P(1)(y’-,,y)exp(t/y  ’  )du’J  Ibu)(t)dt  (3. IS) 

where  the  terms  represent  the  surface  emission,  the  reflection 
of  the  energy  incident  on  the  layer  after  being  attenuated  by 
the  layer  and  the  reflection  of  the  attenuated  gaseous  emission, 
respectively.  Although  Eq.  (3.15)  has  been  formulated  in  gene¬ 
ral,  bidirectional  reflectance  data  is  very  seldom  available  for 
engineering  surfaces  and  the  diffusely  and  specularly  reflecting 
models  are  often  resorted  to.  When  the  surface  is  a  diffuse 
emitter  and  reflector,  eu(p)  ■  eu  and  pu(u'-»u)  ■  pv,  the  into 
gration  over  direction  in  Eq.  (3.15)  can  be  readily  performed 
to  give  the  familiar  exponential  integral  form  (1).  For  a 
specular  reflection,  the  bidirectional  reflectance  is 

Pu(y) 

p^y'-y)  •  2 u « "  6(y’-y) 

where  pw(u)  is  the  directional  reflectivity  and  6(u*-u)  the 
Oirac  delta  function.  The  intensity  leaving  the  plate  then 
becomes 


p«(l,)  '  ' 


,/u 


Ui 


♦ 


I. 


ou 

Pu(y)Ihw(t)exp(-t/u)dt/u. 
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(3.  lb) 


In  order  to  obtain  the  total  radiant  heat  flux,  Liq.  (3.13) 
nas  to  be  integrated  over  the  spectrum.  For  this  integration 
the  total  band  absorptance 

A ( t )  =  [  I  l-cxp(-K  t)ldw  (3.17) 

jAajL  J 

is  introduced.  A  correlation  for  A(t)  based  on  an  exponential 
representation  of  the  mean  line  intensity  and  the  statistical 
band  model  have  been  developed  [34,  35]  for  the  vibrational- 
rotational  bands  of  several  gases.  As  pointed  out  in  the  in¬ 
troduction,  several  researchers  have  used  this  model  and  have 
found  it  to  be  quite  adequate  for  energy  transfer  predictions. 

To  facilitate  the  use  of  the  total  band  absorptance  an 
approximation  for  the  exponential  integral  function  of  the 
form, 

n 

E2(t)  -  l  a.  exp(-b.t)  (3.18) 

i  - 1  1  1 

which  represents  an  approximation  to  the  angular  distribution 
of  the  intensity  is  employed.  Equation  (3.18)  is  a  generaliza¬ 
tion  of  the  exponential  kernel  approximation  [35]  which  can  be 
shown  to  ho  equivalent  to  the  Mi lne-Eddington  approximation  [37]. 
ilowever,  retaining  additional  terms  is  equivalent  to  fourth  and 
higher  order  spherical  harmonics  expansion  of  the  intensity, 
lor  this  study  a  two- term  approximation  with  coefficients 
derived  by  Murty  [37]  was  used  to  achieve  an  improved  acciu-acy. 

Substitution  of  the  exponential  integral  approximation  and 
the  total  band  absorptance  into  Eq.  (3.13),  and  a  subsequent 
integration  over  the  spectrum  leads  to  an  expression  for  the 
total  radiant  flux 

F+(y)  *  ~T  l  {*[  Ioc(0,p)[Aa,-A(y/u)]ydu 
a,r  bands1  K  *c 

*ec^  [A«-A^jjZj]|idii 

+  f  sign(y-z)Ebc(z)J[a^A,(bi|y-z!)dz|.  (3.  19) 

Jo  i  ' 

liquation  (3.19)  was  obtained  assuming  that  the  absorption  bands 
are  narrow  enough  that  the  Planck  function  could  be  considered 
constant  over  the  entire  band  width.  Similarly,  an  express! 
for  the  intensity  emerging  from  the  layer  within  a  given  spe 
tral  band  is 

I£c(«,m)  =  l£c(0,y)!A«-A(6/y)]  ♦  I  Ibc(z)Ia.A* |b. (6-z)/p ^dz 

0  i  (3. 
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3.1.4  Method  of  Solution 


When  coupled  to  radiation  transfer  the  governing  equations 
represent  a  system  of  two-point,  nonlinear,  integro-di fferenti al 
equations  for  which  a  closed-form  solution  is  not  possible  and 
numerical  methods  must  be  resorted  to.  Numerically,  coupled 
integral  equations  are  simpler  to  solve  than  coupled  integro- 
differential  equations.  An  integration  of  the  non-dimensional 
conservation  relations  reduces  them  to  a  system  of  integral 
equations 

P*  Cx)  '  1  *  f  *  jfr  (8„-S)]dX  (3.21) 

»,(x)  -  *  f  U-w.)  -  3|] dx  (3.22) 

0(X)  =  Q00+J  Bo  F++  j+ACp  +  NM+j+  0  dX  (3.23) 

ok  v  *  r  ' 


where  x  *  1  -  £.  The  numerical  solution  was  accomplished  by 
successive  approximations  with  a  weighted  averaging  of  the 
previous  solutions.  It  was  found  that  in  the  first  approxima¬ 
tions  the  weighting  of  the  current  approximation  should  be 
small  to  assure  reasonably  rapid  convergence. 


3.2  Experiments 


3.2,1  Description  of  Experimental  Facility 


The  test  apparatus  was  installed  in  a  fume  hood  in  order 
to  guarantee  a  sufficient  degree  of  environmental  control.  A 
fan  drew  air  into  the  hood  from  below  the  chamber  and  exhausted 
it  through  a  duct  into  the  atmosphere.  A  slightly  reduced 
pressure  in  the  fume  hood  prevented  the  transpired  gases  from 
endangering  the  personnel  and  also  provided  a  means  of  control¬ 
ling  the  flow  field  of  the  injected  gas.  The  test  apparatus 
consisted  of  the  traversing  mechanism,  support  stand,  and  test 
section.  A  photograph  of  the  test  apparatus  is  shown  in  Fig.  9. 


The  traversing  mechanism  was  designed  and  constructed  by 
Research  Services  Inc.,  St.  Paul,  Minn.  The  radiation  detector 
was  mounted  on  a  yoke  shaped  scanning  mechanism  permitting  the 
positioning  of  the  detector  at  various  polar  angles.  The  arms 
of  the  yoke  were  adjustable  such  that  the  distance  between  the 
radiation  detector  and  the  axis  of  rotation  could  be  varied 
from  4"  to  24"  at  spacing  intervals  of  one  inch.  The  detector 
holder  could  be  displaced  along  the  yoke  (6"  to  the  left  of  the 
center  line  and  6"  to  the  right).  An  angle  indicator  provided 
with  locking  pins  served  for  variation  of  polar  angles  (from 
the  vertical)  in  S°  increments.  The  rotating  component  of  the 
scanning  mechanism  was  mounted  on  a  frame  which  could  be  moved 
longitudinally  by  about  12"  sliding  on  rectangular  rods. 

Besides  serving  for  the  positioning  of  the  detector  the 
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traversing  mechanism  was  also  utilized  for  the  support  of  the 
aspirating- temperature  probe.  The  test  section  was  installed 
in  the  traversing  mechanism  such  that  its  axis  of  rotation 
always  remained  in  the  surface  of  the  porous  plate,  and  the 
longitudinal  and  lateral  motions  were  always  parallel  to  the 
plate  surface. 

3.2.2  Test  Assembly 

The  heated  porous  plate  was  housed  in  a  12"  x  9"  x  4" 
sheet  metal  box  open  at  one  end.  Welded  to  the  chamber  were 
angle  irons  which  served  to  attach  it  to  the  support  frame. 

The  test  section  consisted  of  a  gas  plenum  chamber  closed  on 
one  side  by  the  porous  plate,  copper  electrodes,  transite 
insulating  material,  and  sheet-metal  gas  layer  extenders.  The 
test  section  is  shown  in  Fig.  10. 

A  porous  plate  with  large  internal  surface  area  and  uni¬ 
form  porosity  is  desired  for  effective  heat  transfer  and 
flow  distribution.  A  survey  of  available  materials  revealed 
that  a  nickel-base  permeable  alloy  sheet,  called  Poroloy,  pro¬ 
duced  by  the  Bendix  Filter  Division  of  the  Bendix  Corporation, 
was  appropriate.  The  Poroloy  sheet  was  composed  of  a  nickel 
base  with  224  chromium,  18.5%  iron,  9%  molybdenum,  1.5%  cobalt, 

6%  tungsten,  and  small  percentages  of  carbon,  manganese,  silicon, 
phosphorus,  and  sulfur.  This  alloy  is  oxidation  resistant  up  to 
2100“F  and  has  a  relatively  high  strength  above  1500°1;.  The 
Poroloy  plates  were  made  from  a  0.005  inch  diameter  wire, 
flattened  to  a  ribbon,  wound  into  a  laminated  structure  at  a 
nominal  30  degree  wind  angle  (60  degree  crossover  angle),  and 
sintered. 

Two  6"  x  12"  sheets  of  different  permeabilities  were 
acquired.  One  sheet  has  a  permeability  of  2.0  x  10'®  ft  and 
the  other  with  a  permeability  of  2.0  x  10'®  ft  was  ten  times 
more  porous  than  the  first.  The  permeability  of  the  individual 
sheets  was  tested  by  the  manufacturer  to  conform  to  within  ±  10% 
of  the  level  specified  when  subjected  to  a  pressure  difference 
of  10.0  psi  with  exit  pressure  equal  to  14.7  psia  and  air 
temperature  of  70°F. 

Due  to  the  alloy’s  relatively  low  electrical  resistance 
(approximately  6.9  x  10'1  ohm  at  70°F  for  a  2.2"  x  3.5"  x  .025" 
sheet),  a  current  of  ISO  to  250  amperes  supplied  by  a  DC  arc 
welder  was  needed  to  reach  the  desired  plate  temperature.  In 
order  to  minimize  heat  conduction  losses  from  the  plate  and 
still  maintain  structural  strength,  transite  pieces  were  placed 
in  the  box. 

The  upstream  part  of  the  sheet-metal  box  consisted  of  a 
plenum  chamber.  A  gas  line  from  a  battery  of  supply  bottles 
and  various  measuring  and  regulation  devi ccs  were  connected  to 
the  top  of  the  plenum  with  the  gas  flowing  into  the  chamber 
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Figure  10.  Sectional  View  of  Test  Assembly 
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and  out  through  the  porous  plate.  A  deflector  plate  was 
mounted  between  the  porous  plate  and  the  gas  inlet  in  order 
to  guarantee  uniform  gas  flow  through  the  plate  and  to  shield 
the  plenum  from  radiant  heating.  Attached  to  the  deflector 
plate  was  a  flexible  steel  conduit  which  served  as  an  instru¬ 
mentation  port,  allowing  surface  and  gas  temperature  measure¬ 
ments  inside  the  plenum. 

Sheet  metal  guards  of  varying  length  could  be  mounted 
below  the  plate  to  serve  as  gas  layer  extenders.  The  gas 
layer  extenders  were  thermally  insulated  on  both  sides  with 
asbestos  material  in  order  to  minimize  heat  conduction  losses. 
The  direction  of  the  gas  flow  was  downward  until  it  reached 
the  exit  plane  of  the  extenders  where  buoyancy  reversed  its 
flow  direction. 

3.2.3  Instrumentation  and  Calibration 

In  order  to  predict  the  radiation  emitted  by  the  plate 
its  surface  temperature  and  the  local  variation  as  well  as  the 
emittance  have  to  be  known.  A  Leeds  and  Northrup  surface 
temperature  probe  ("Surtemp")  was  used  to  measure  the  tempera¬ 
ture  distribution.  Its  specifications  were: 

Accuracy:  Ambient  to  530°F:  ±  3°F;  530°  to  1200°F: 

±  (3/81  of  reading  +  1°)F. 

The  "Surtemp"  probe  had  an  absolute  error  of  5°F  at  1000°F. 
Since  the  probe  had  to  be  operated  manually  (i.e.,  it  had  to 
be  held  in  position  touching  the  lower  surface  of  the  plate) 
and  could  not  be  used  while  radiation  measurements  were  being 
taken,  chromel-alumel  surface  thermocouples  were  employed 
during  certain  phases  of  the  experiments. 

Previous  experience  accumulated  by  research  personnel  at 
the  Allison  Division  of  the  General  Motors  Corporation  was 
utilized  for  the  design  and  construction  of  an  aspirating  gas 
temperature  probe  for  the  measurement  of  gas  layer  temperature. 
The  probe  consisted  of  Inconel  tubing  with  an  insulated  0.02" 
O.D.  Inconel  tube  containing  the  thermocouple  wires,  see  Fig. 
11.  The  thermocouple  wires  were  of  36  gage  (0.003"D)  chromel- 
alumel.  A  tee  with  two  air-line  connections  was  positioned  at 
the  base  of  the  tubing  with  one  opening  smaller  than  the  other. 
As  air  was  blown  into  the  small  opening  and  out  of  the  large 
opening,  the  resulting  jet  lowered  the  pressure  of  the  tube 
below  atmospheric  pressure  causing  the  surrounding  gas  to  be 
drawn  into  the  probe  through  a  5/64"  D  hole  in  the  tip.  The 
tip  was  normal  to  the  flow  so  that  the  hole  faced  in  the  flow 
direction.  The  thermocouple  junction  was  made  by  a  small  head 
weld  (0.010"0)  and  was  shielded  by  a  Monel  tube. 

The  transpiring  gas  flow  was  measured  with  a  Fisher-Porter 
Rotameter  which  was  calibrated  against  a  standard  integrated 
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Figure  11.  Sectional  View  of  Gas  Temperature  Probe  Tip 


gas  flow  meter  manufactured  by  American  Meter  Company.  The 
two  flow  meters  were  mounted  next  to  each  other  and  connected 
in  such  a  way  that  the  flow  through  the  rotameter  could  be 
directed  through  the  integrating  flow  meter.  Once  the  rotameter 
was  calibrated,  the  integrating  flow  meter  was  by-passed  to 
avoid  an  unnecessary  pressure  drop. 

Power  dissipation  was  determined  from  the  product  of  the 
input  voltage  and  current.  The  voltage  drop  across  the  test 
section  was  measured  with  a  standard  voltmeter.  Electric 
current  was  determined  by  measuring  the  voltage  drop  across  a 
shunt  resistance  of  125  w  ohms. 

The  radiant  energy  emerging  from  the  gas  layer  was  measured 
with  an  infrared  radiometer  which  consisted  of  a  wide-band-pass 
filter,  a  potassium-bromide  lens,  two  diaphragms,  and  a  detector, 
(see  fig.  12).  The  detector,  a  thermopile,  was  manufactured  by 
Charles  M.  Reeder  Company  of  Detroit,  Michigan.  The  filter  and 
1/4”  l)  diaphragm  were  mounted  in  front  of  the  thermopile.  A 
second  3/8"  P  diaphragm  was  positioned  in  front  of  the  lens  and 
both  were  held  in  a  1/2"  D  tube.  The  thermopile  was  then 
inserted  into  the  tube  behind  the  lens.  The  assembled  radiometer 
was  mounted  and  adjusted  on  the  traversing  mechanism  so  that  any 
desired  surface  clement  of  the  plate  could  be  brought  into  focus. 
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In  order  to  measure  the  attenuation  in  a  particular  wave¬ 
length  band  a  filter  was  used.  For  example,  the  filter  employed 
for  the  4.3y  (fundamental  band  of  C02)  had  a  band  center  at 
4.35y  and  a  half  band  width  of  0.43y.  A  transmittance  curve  is 
shown  in  Fig.  13.  The  filter,  0.040"  thick,  was  manufactured 
by  the  Optical  Coating  Laboratory,  Santa  Rosa,  California. 

In  order  to  investigate  the  attenuation  due  solely  to  the 
presence  of  participating  gas,  the  emerging  radiant  energy  was 
compared  with  that  measured  with  a  transparent  gas  (air)  present. 
Thus  the  effect  of  transpiration  cooling  on  the  emitted  radiation 
was  practically  eliminated.  The  measured  normalized  emergent 
band  intensity  was  expressed  as 

‘exp'*-8)  *  lx  V4-8)dl/j1A.  *»  'x'0'8)-11  (5-24» 

This  definition  is  analogous  to  that  of  Eq.  (2.S),  except  that 
the  spectral  transmittance  of  the  optical  system  t^  (wide  band 
pass  filter  together  with  the  KBr  lens  and  thermopile  window, 
see  Fig.  12)  has  been  incorporated. 

3.2.4  Calibration  and  Test  Procedure 

In  order  to  determine  the  effective  responsivity  of  the 
radiometer,  it  was  necessary  to  evaluate  the  total  transmittance 
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of  the  filter,  the  focal  length  of  the  lens,  and  the 
wavelength  band  responsivity  of  the  thermopile.  The 
of  the  total  radiant  energy  coming  from  a  black  body 
by  a  small  solid  angle  and  passing  through  a  filter 


total  and 
fraction 
subtended 
is  given  by 


f  tFXIbx(T)dX/(aT4/iO 


(3.25) 


The  results  were  computed  and  found  to  be  0.0573  for  T  =  1000°F, 
0.0560  for  T  =  1250°F,  and  0.0514  for  T  =  1500°F.  The  focal 
length  of  the  lens  was  determined  by  two  different  methods  [38]; 
both  methods  lead  to  a  focal  length  of  0.866. 

In  order  to  determine  the  responsivity  of  the  thermopile, 
a  calibration  curve  was  experimentally  established.  Since  each 
thermopile  had  different  characteristics  it  is  important  that 
this  measurement  be  done  with  extreme  care  for  each  thermopile. 
The  detector  was  irradiated  with  a  calibrated  standard  black  body 
source  (Infrared  Industries  Model  No.  406)  as  recommended  by 
Nicodemus  [39] .  The  output  signal  of  the  thermopile  was  recorded 
as  a  function  of  incident  radiant  energy.  The  black  body  temper¬ 
ature  was  controlled  with  a  calibrated  temperature  controller. 

The  control  range  extended  from  50°C  to  1000°C  with  a  specified 
accuracy  of  ±  1°C. 

In  essence  two  types  of  calibration  curves  were  determined 
for  the  radiometer.  The  first  type  pertains  to  the  thermopile 
alone  on  a  total  energy  basis  and  the  second  to  the  thermopile 
together  with  the  filter  on  a  band  basis.  Both  types  of  cali¬ 
bration  curves  were  established  for  source  temperatures  of 
1000°F,  1250°F,  1500°F.  One  of  these  curves  is  illustrated  in 
Fig.  14  and  shows  that  the  responsivity  of  the  thermopile  on  a 
total  energy  basis  is  linear. 

To  determine  the  effective  transmittance  of  the  radiometer 
without  a  filter  the  lens  and  its  tubular  retainer  were  mounted 
at  a  fixed  distance,  s,  from  the  blackbody  source.  The  image 
plane  of  the  lens  was  located  by  adjusting  the  distance  between 
the  thermopile  and  lens  for  the  maximum  thermopile  output.  From 
the  black  body  temperature,  geometry  and  thermopile  calibration, 
the  effective  transmittance  could  be  calculated.  Transmittances 
as  a  function  of  the  two  independent  geometric  variables,  object 
distance,  s,  and  black  body  source  diameter  (the  remainder  of 
the  geometry  is  given  by  the  thin  lens  formulas)  were  obtained. 

Voltage  signals  from  the  thermocouples  and  thermopile  were 
connected  by  chromel-alumel  couplings  and  leads  to  an  ice  bath 
reference  junction.  The  signal  was  recorded  through  a  10  channel 
switching  mechanism  on  either  a  Leeds  and  Northrup  type  K-3 
potentiometer  or  s  Honeywell  Visicorder  data  acquistion  system. 

After  the  equipment  was  properly  arranged  and  aligned,  the 
arc  welder  was  ignited  and  the  current  adjusted  to  the  desired 
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Figure  14.  Thermopile  Calibration  Curve 


level.  The  plate  was  allowed  to  reach  the  steady  state,  which 
required  between  one  and  one  half  to  two  hours  of  time.  Sub¬ 
sequently,  the  gas  flow  was  started  and  adjusted  to  the  desired 
rate  and  the  fan  damper  on  the  fume  hood  set.  Final  adjustments 
were  made  in  the  current  level  to  reach  the  required  temperature 
of  the  plate. 

After  the  system  reached  steady  state  conditions  the  volt¬ 
ages  of  the  thermocouples  mounted  on  the  plate  surface  were 
monitored  and  recorded.  Simultaneously  the  plate  surface  tem¬ 
perature  distribution  was  measured  with  the  "Surtemp"  probe. 

The  temperature  of  the  gas,  its  flow  rate,  and  test  section 
pressure  were  also  recorded.  The  gas  temperature  probe  was 
positioned  1”  below  the  gas  layer  extender  while  the  air  flow 
through  the  expansion  valve  was  recorded.  To  verify  the  steady 
state  the  data  were  monitored  with  the  Visicorder  with  the 
final  data  being  obtained  with  the  K-3  potentiometer. 
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Radiation  measurements  were  also  made  under  steady  state 
conditions.  Since  the  temperature  distribution  of  the  plate 
changed  with  flow  rate,  it  became  necessary  to  compare  the 
emerging  radiation  of  an  essentially  transparent  gas  (air)  to 
that  of  an  absorbing-emitting  gas  under  otherwise  identical 
conditions.  The  injected  air  flow  rate  was  set  and  the  tempera¬ 
ture  distribution  of  the  plate  measured  with  the  "Surtemp"- 
probe.  Once  radiation  data  were  obtained  with  air  as  the  work¬ 
ing  fluid,  it  was  replaced  by  a  radiation  participating  gas 
udner  corresponding  conditions. 

The  signal  from  the  radiometer  was  strong  enough  so  that 
it  could  be  read  with  a  potentiometer  without  further  amplifi¬ 
cation.  During  testing  it  was  found  that  the  temperature  of 
the  thermopile  casing  rose.  In  order  to  compensate  for  the 
drift  of  the  thermopile  output,  a  metal  plate  of  uniform  tem¬ 
perature  was  placed  between  the  source  and  the  radiometer  and 
the  output  of  the  thermopile  measured.  The  plate  was  subse¬ 
quently  removed  and  the  signal  corresponding  to  the  radiation 
received  by  the  thermopile  was  measured  as  quickly  as  possible. 
Subsequently  the  metal  shield  was  returned  to  the  radiation 
path  in  order  to  keep  the  radiometer  at  room  temperature. 

3 • 3  Results  and  Discussion 

3.3.1  Analytical  Results 

Since  the  independent  parameters  are  numerous,  a  review 
of  them  seems  appropriate.  They  are  the  Re,  Gr,  Sc,  Pe,  N, 
and  Bo  numbers,  the  gas  temperature  next  to  the  plate,  the 
radiation  characteristics  of  the  plate,  the  gas  absorption 
properties,  the  incident  radiation,  and  the  optical  or  layer 
thickness.  Results  are  presented  for  the  diffusely  emitting 
and  reflecting  surface  in  the  absence  of  external  radiation  in 
order  to  reduce  the  number  of  parameters.  Carbon  dioxide, 
methane,  and  sulfur  dioxide  were  selected  for  the  computations 
because  of  their  strong  fundamental  bands  in  the  near  infrared. 

The  results  are  presented  in  terms  of  a  normalized  emergent 
band  intensity  defined  by,  see  Eqs.  (2.4)  and  (2.S), 

I*  *  f  Iu(«,u)d«/|  lwt(0,y)dw  (3.26) 

where  lut(0,u)  is  the  intensity  leaving  the  plate  in  the  absence 
of  an  absorbing  gas  layer.  The  normalized  band  intensity  1* 
serves  as  a  measure  of  the  attenuation.  A  certain  amount  of 
ambiguity  is  associated  with  an  effective  band  width  for  an 
exponential  type  band  model.  At  high  pressure,  the  exponential 
wide  band  model  is  (34) 

*u(T)  -  pC*  exp(-|w  -  wc|/A0)  (3.27) 
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where  A„  is  a  bandwidth  parameter  which  depends  directly  on 
temperature  and  C\  is  the  line  intensity  to  spacing  ratio  which 
varies  inversely  with  temperature  for  typical  gas  bands.  A0  , 
the  hand  half-width  could  be  used  as  a  representative  bandwidth, 
hut  would  he  too  narrow  to  be  meaningful.  From  Eq.  (3.27)  an 
effective  bandwidth  may  be  defined  as 

Aw  =  2A0  In  [<(wco)/<(“c)j  (3.28) 

where  wc  is  the  band  cut-off  wave  number.  The  band  cut-off  to 
hand  center  absorption  coefficient  ratio  k(wco)/k(wc)  was  taken 
as  0.001  in  order  to  remain  consistent  with  the  carbon  monoxide 
effective  bandwidth  presented  by  Tien  [3]. 

The  effect  which  the  blowing  rate  has  on  the  temperature 
and  concentration  profiles  is  shown  in  Figs.  15  end  16.  A 
measure  of  the  T-D  and  D-T  effects  and  the  effect  of  coupled 
radiative  transfer  can  be  obtained  from  a  comparison  with  an 
analysis  neglecting  these  effects.  From  the  concentration 
profile,  it  is  seen  that  the  T-D  and  D-T  effects  are  quite 
weak  and  might  well  be  neglected.  The  coupling  of  the  radiative 
transfer  is  also  quite  small  because  of  the  low  partial  pressures 
of  the  absorbing  gas  and  the  fact  that  only  about  3%  of  the  radi¬ 
ant  energy  emitted  by  the  plate  lies  within  the  absorbing  bands. 
It  is  of  interest  to  note  that  even  with  the  weak  coupling,  the 
trends  in  the  temperature  profiles  are  the  same  as  those  for 
stronger  coupling.  Both  the  temperature  and  concentration  in¬ 
crease  with  increasing  blowing  rate.  Since  the  concentration 
is  increasing,  the  gas  layer  is  becoming  more  opaque  within  the 
absorption  bands  and  the  energy  leaving  the  plate  is  more 
effectively  attenuated.  The  net  effect  is  one  of  reducing  the 
emergent  intensity  as  shown  in  Fig.  17. 

A  second  effect  is  due  to  the  variation  of  the  total  band 
absorptance  parameters.  Neglecting  the  pressure  broadening, 
the  total  band  absorptance  is  a  function  only  of  the  integrated 
band  intensity  and  bandwidth.  For  carbon-dioxide,  the  bandwidth 
variation  is  about  the  same  for  both  bands,  while  the  integrated 
band  intensity  is  one  order  of  magnitude  higher  for  the  funda¬ 
mental  4.3u  band.  The  trends  of  the  15y  band  which  is  a  rela¬ 
tively  broad  band  are  primarily  dependent  on  the  bandwidth, 
whereas  the  4.3p  band  is  primarily  dependent  on  the  integrated 
band  intensity.  Both  parameters  determine  the  attenuation  of 
radiation  by  the  gas  layer,  but  from  the  results  in  Fig.  17  tli 
integrated  band  intensity  is  the  more  dominant.  This  effect 
was  also  observed  in  calculations  in  which  the  gas  temperature 
next  to  the  plate  was  varied  while  the  plate  temperature  was 
maintained  constant. 

The  attenuation  as  a  function  of  the  plate  emittance  with 
polar  angle  and  bandeenter  wavelength  as  parameters  is  present  d 
in  Fig.  18.  Although  the  temperature  and  concentration  distri¬ 
butions  are  not  shown,  they  were  computed  and  found  to  vary  only 
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Figure  15.  Effect  of  Injection  Rate  on  the  Temperature 

Distribution;  C02 ,  6  =  6cm,  Tp  =  836°K,  e  =  0.6 
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Effect  of  the  Plate  Emittance  on  the  Normalized 
Emergent  BanU_ Intensity ;  C02,  Re  =  2,  6  =  2cm, 


slightly  with  plate  emittance.  Inspection  of  Fig.  13  reveals 
that  there  is  a  plate  emittance  at  which  the  actually  emergent 
intensity  in  the  presence  of  an  absorbing  gas  equals  that  in 
presence  of  a  non- absorbing  gas.  This  behavior  will  be  referred 
to  as  the  "break-even  emittance."  As  the  plate  emittance  is 
decreased,  the  gas  emission  which  emerges  directly  from  the 
layer  remains  almost  constant  and  the  energy  reflected  by  the 
plate  increases.  It  follows  that  the  emergent  intensity  of  a 
surface  with  an  absorbing  gas  present  can  exceed  that  in  absence 
of  an  absorbing  gas  at  low  plate  emittance  values.  An  attenua¬ 
tion  of  the  energy  emitted  from  the  plate  takes  place  only  when 
the  emittance  is  greater  than  the  "break-even  emittance." 

The  dependence  of  the  emergent  intensity  on  the  plate  tem¬ 
perature  and  layer  thickness  is  presented  in  Fig.  19.  As 
expected,  an  increase  of  the  layer  thickness  enhances  the  atten¬ 
uation.  Inspection  of  Fig.  19  reveals  further  that  the  attenua¬ 
tion  is  approaching  a  limit.  In  the  4.3p  band  the  results  for 
a  12  cm.  layer  thickness  were  within  one  plotted  linewidth  of 
those  for  an  8  cm.  layer.  At  this  limit,  the  gas  becomes  opaque 
in  the  band  center.  Energy  emitted  by  the  plate  continues  to 
emerge  from  the  layer  through  the  band  wings  which  do  not  become 
opaque  as  rapidly  as  the  band  center.  The  limit  is  thus  the 
layer  thickness  at  which  the  attenuation  by  the  band  center  has 
reached  an  optimum  and  any  further  increases  in  the  attenuation 
must  be  accomplished  by  the  band  wings. 

In  the  near  normal  direction,  the  emergent  intensity  is 
almost  isotropic  since  the  path  length  does  not  change  signifi¬ 
cantly.  At  oblique  angles,  the  attenuation  increases  quite 
rapidly  because  of  the  increase  in  path  length.  Results  were 
not  obtained  for  0  >  80°  because  of  the  singularity  in  Eq. 

(3.22)  as  u  =  cos0  approaches  zero.  As  0  -*•  it/2,  I*  -*■  0  which 
is  consistent  with  the  results  presented  in  Figs.  19  and  20. 

To  investigate  the  influence  of  the  absorption  properties 
of  the  injected  gas,  calculations  using  the  absorption  charac¬ 
teristics  of  methane  [40]  and  sulfur  dioxide  [41],  considering 
only  their  strongest  absorption  bands,  were  made.  To  make  the 
results  compatible  with  those  obtained  for  carbon-dioxide,  the 
Reynolds  number  was  adjusted  accordingly.  The  results  for 
methane  and  sulfur  dioxide  with  the  layer  thickness  as  a  param¬ 
eter  are  presented  in  Fig.  20.  Since  sulfur  dioxide  is  the 
stronger  absorber  in  the  7.35p  fundamental  band  than  the  4 . 3u 
fundamental  CO2  band,  the  attenuation  increases,  especially  for 
the  smaller  layer  thicknesses.  Hence,  SO2  would  be  advantageous 
in  physical  systems  with  small  boundary  layers.  The  fundamental 
SO 2  and  C1U  bands  are  also  wider  than  the  fundamental  C02  hands 
and  therefore  cover  a  larger  portion  of  the  spectrum.  Although 
SO 2  is  a  stronger  absorber  and  attenuates  a  larger  portion  of 
the  spectrum,  its  fundamental  band  is  remote  from  the  maximum 
plate  emission  for  the  temperatures  considered.  The  net  amount 
of  energy  attenuated  by  the  SO2  fundamental  band  is  consequently 
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not  as  great  as  that  attenuated  by  the  fundamental  band  of  C02. 
On  a  total  basis,  the  4 . 3u  CIU  band  could  be  more  effective 
than  the  SO?  fundamental  band. 

3.3.2  Experimental  Results 

The  experimentally  measured  centerline  gas  temperature 
profiles  with  C02  gas  injected  through  the  porous  plate  are 
shown  in  Fig.  21  with  blowing  Reynolds  number  as  a  parameter. 

The  temperature  profiles  differ  appreciably  from  those  that 
would  be  expected  for  one-dimensional  flow  and  heat  transfer, 
particularly  near  the  exit  plane  of  the  gas  layer  extenders. 

This  is  partly  due  to  effective  mixing  of  the  gas  streams.  It 
was  observed  that  the  gas  temperature  in  the  exit  plane  could 
be  controlled  by  the  amount  of  air  drawn  through  the  hood.  The 
air  flow  in  the  hood  was  so  adjusted  that  the  effect  of  the 
mixing  region  on  the  temperature  distribution  in  the  gas  layer 
was  as  small  as  possible. 

The  radiation  attenuation  results  are  given  in  Tables  3, 

4  and  5  and  Fig.  22.  The  normal  spectral  emissivity  of  the 
plate  increased  with  the  level  of  oxidation.  The  emissivity 
of  the  Poroloy  plate  (see  Fig.  23)  was  measured  [42]  and  found 
to  be  a  function  of  temperature  and  oxidation  level.  A  well 
oxidized  plate  was  used  with  a  corresponding  spectral  normal 
emissivity  of  approximately  0.65  in  the  wavelength  band  of  4.3y. 

The  effect  of  the  blowing  Reynolds  number  on  attenuation 
was  found  to  be  small.  The  trend  indicates  an  increase  in 
attenuation  with  increasing  blowing  Reynolds  number,  but  the 
data  were  not  consistent  enough  to  be  conclusive.  Increasing 
layer  thickness  and  plate  temperature  enhanced  the  attenuation. 
The  effect  of  the  viewing  angle  was  slight  since  only  near 
normal  directions  could  be  observed  (9  <  20°).  The  polar  angle 
affects  the  attenuation  in  two  ways.  First,  the  emissivity  of 
the  plate  is  a  function  of  direction;  however,  the  directional 
emissivity  of  metals  does  not  vary  appreciably  with  direction 
angles  less  than  45°  [10].  Second,  the  optical  path  length  of 
the  radiation  propagating  through  the  gas  is  inversely  propor¬ 
tional  to  the  cosine  of  the  angle  of  observation,  see  Eq.  (2.1). 
For  9  *  20°,  l/cos8  ■  1.06,  a  value  which  indicates  that  the 
path  length  has  not  increased  appreciably  and  hence  no  signifi¬ 
cant  change  is  expected. 

3.3.3  Comparison  of  Analysis  and  Experiment 

The  predicted  temperature  profiles  differed  appreciably 
from  the  measured  ones.  This  discrepancy  was  especially  large 
in  the  vicinity  of  the  exit  of  the  gas  layer  extenders.  The 
reasons  for  this  difference  are  clear.  In  the  analysis  the 
temperature  and  concentration  at  the  edge  of  the  layer  (y  ■  6) 
were  imposed,  and  the  mixing  region  of  the  two  opposing  streams 
was  assumed  to  be  negligibly  thin;  however,  flow  visualization 
revealed  extensive  mixing  in  this  region  and  that  it  extended 
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curve  Re  8 ,  cm. 

1  43.5  7.62 

2  29  5.06 

3  130.5  7.62 
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Table  3 .  Effect  of  Blowing  Reynolds  Number  on  Attenuation 


for  6  =  7.62  cm,  T  =  822°K,  and  6=0° 


Rc 
I  * 

exp 


4  3.5 
0.502 


97.0  130.5 

0.521  0.512 


Table  4.  Effects  of  Plate  Temperature  and  Gas  Layer  Thickness 
on  Emergent  Band  Intensity  for  9=0° 


6  ( cm) 

V  °K 

822 

960 

1090 

5.08 

0.437 

0.481 

0.480 

7.62 

0.471 

0.525 

0.502 

10.16 

0.492 

0.557 

0.542 

Table  5.  Effect 

6  ■  5 . 

of  Viewing 

08  cm,  T  ■ 
P 

Angle  on  Attenuation 

960°K,  and  Re  «  29 

for 

9(°) 

20 

0 

-20 

'exp 

0.482 

0.500 

0.490 
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liffect  of  Gas  Layer  Thickness  on  Attenuation 
in  ■  0.01  lb/min,  6  *  0 


x  0 


below  the  exit  plane  of  the  gas  layer  curtain.  Also  in  the 
analysis  lateral  heat  conduction  was  neglected.  Although  con¬ 
duction  losses  to  the  sides  were  reduced  by  insulating  the 
metallic  extenders,  some  energy  transfer  in  the  lateral  direction 
could  not  be  avoided. 

The  measured  attenuation  is  larger  than  the  predicted  one. 
There  are  a  number  of  reasons  for  this  difference.  The  extension 
of  C02  concentration  beyond  the  exit  plane  of  the  gas  layer  cur¬ 
tain,  as  already  discussed,  would  increase  the  optical  path  and 
hence  the  attenuation.  In  addition,  the  filter  band  width  used 
with  the  radiation  detector  was  somewhat  smaller  than  the  effect¬ 
ive  band  width  defined  in  the  analysis,  Eq.  (3.28),  and  hence  a 
smaller  attenuation  has  to  be  expected  in  the  predictions.  Also 
previous  investigators  have  shown  that  the  exponential  wide  band 
model  for  the  total  band  absorptance  underpredicts  the  total 
band  absorptance  for  nonisothermal  systems  [43].  The  introduc¬ 
tion  of  an  effective  bandwidth  was  necessitated  because  the 
spectral  absorption  coefficient  is  a  function  of  temperature, 
see  Eq.  (3.27).  It  should  be  emphasized  here  that  the  differ¬ 
ence  in  the  two  definitions  of  the  attenuation  is  caused  by  the 
fixed  bandwidth  of  the  filter  as  opposed  to  varying  bandwidth  in 
the  analytical  model.  There  was  no  possible  logical  and  consist¬ 
ent  way  to  adjust  the  correlation  constants  and  the  bandwidth  in 
the  total  band  absorptance  equation  [3,  34,  35]  to  correspond 
with  the  bandwidth  of  the  wide  band  pass  filter  used  with  the 
radiation  detector.  For  these  reasons  a  direct  comparison  of 
experimental  and  analytical  results  remains  somewhat  ambiguous; 
however,  the  trends  are  correctly  predicted  by  the  analytical 
model . 
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4.  ATTENUATION  OF  RADIATION  BY  AN  ABSORBING 
GAS  FLOWING  OVER  A  HOT  PLATE 


4 .  I  Analysis 

4.1.1  Physical  Model  and  Assumptions 

The  transpiration  cooling  of  critical  structural  components 
of  jet  engines,  rocket  motors,  airfoil  surfaces,  etc.,  which  are 
subjected  to  high  temperatures,  is  an  effective  cooling  technique 
and  has  received  considerable  attention  [9],  This  can  be  done 
by  designing  the  porosity  of  the  surfaces  so  as  to  obtain  a 
desired  distribution,  along  the  walls,  of  normal  mass  flow  of  the 
coolant.  In  addition  to  the  cooling,  if  the  injected  gas  is  also 
capable  of  absorbing  radiation,  the  radiation  leaving  the  hot 
surface  will  be  partly  attenuated  and  would  decrease  the  amount 
of  radiant  energy  incident  on  a  detector  some  distance  away.  The 
main  purpose  of  the  present  analysis  is  to  predict  the  flow  vari¬ 
ables  and  the  attenuation  for  an  injected  gas  which  is  capable 
of  absorbing  and  emitting  radiation  only  in  certain  wavelength 
bands . 

The  physical  model  considered  in  the  analysis  is  shown  in 
Fig.  24.  The  model  shown  in  Fig.  24a  has  a  solid  starting  length, 
Xq ,  followed  by  a  heated  porous  section.  The  solid  plate  up¬ 
stream  of  the  heated  porous  one  is  adiabatic.  The  coolant  gas 
under  uniform  flow  conditions  is  injected  into  the  free  stream 
air  from  the  bottom  of  the  plate.  In  the  model  illustrated  in 
Fig.  24b  the  coolant  is  injected  through  a  cold  porous  plate 
followed  by  a  solid  heated  section  having  a  starting  length  x0. 

The  porous  plate  is  adiabatic  and  the  coolant  is  injected  uni¬ 
formly  into  the  main  stream.  In  both  models  the  porous  and  the 
solid  plates  are  heated  uniformly.  The  following  assumptions 
arc  made  in  the  analysis: 

1.  The  flow  is  laminar  and  two-dimensional  with  zero 
pressure  gradient. 

2.  There  are  no  chemical  reactions  in  the  fluid. 

3.  The  thermal  diffusion  and  diffusion-thermo  effects 
are  negligible. 

4.  The  transport  properties  are  assumed  to  have  a  simple 
power  law  temperature  dependence  and  the  Prandtl  and 
Schmidt  numbers  are  taken  to  be  constant. 
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Figure  24.  Physical  Models  of  the  Problem 
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5.  The  variation  of  the  specific  heat  of  the  coolant 
and  mixture  with  temperature  is  neglected. 

().  The  velocities  are  low  enough  so  that  viscous  heat 
dissipation  effects  can  be  neglected. 

7.  The  analysis  presented  is  restricted  to  porous  walls 
of  sufficient  compactness  to  insure  that  the  fluid 
and  material  temperatures  are  essentially  the  same 
throughout  the  wall  and  equal  at  the  surface. 

8.  The  conduction  in  the  porous  plate  is  one  dimensional. 

The  multicomponent  gas  system  is  treated  as  an  effective 
binary  gas  system  composed  of  the  main  stream  gas  "1"  and  the 
injected  gas  "2".  The  properties  of  the  gas  mixture  within  the 
boundary  layer  are  functions  of  the  local  temperature  and  con¬ 
centration  . 


4.1.2  Conservation  Equations 

The  mathematical  model  for  the  assumed  physical  problem  is 
prescribed  by  the  conservation  equations  of  mass,  momentum, 
energy  and  species.  The  governing  boundary  layer  equations  for 
the  flow  of  a  binary,  radiating  gas  over  a  porous  flat  plate  can 
be  written  as  [1,  10,  32] 


Mass:  -  0 

Momentum:  p(u  57  *  v  3y)  *  If) 

«■  .  J  3h.  dh]  3f.  3T)  3F  A  n  3w 

bnerg>  .  p|u  ♦  v  yy-j  --^(k  ^  ♦  pDu  ^ 

Species  (injected):  p(u  *  v  I7)  *  3y[pD 


Uhrh») 

3y 

3w 

w 


(4.1) 

(4.2) 

(4.3) 

(4.4) 


In  the  formulation  of  Eq.  (4.3),  the  radiative  transfer  in  the 
x-direction  has  been  neglected,  following  the  rationale  of 
Refs.  1  and  10.  The  total  (integrated  over  the  spectrum)  radi¬ 
ative  flux  is  given  in  terms  of  the  monochromatic  flux,  Eq. 
(3.13),  through  the  implicit  definition  F  •  /#  F^dX. 

The  boundary  conditions  imposed  are  as  follows: 

Momentum:  u  «  0,  v  *  vw  (a  constant)  at  y  ■  0 

u  -*■  ub ,  v  0  as  y  ♦  *  (4.5) 
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In  transpiration  cooling  the  local  heat  transfer  and  the 
rate  of  injection  of  the  coolant  are  not  independent  but  are 
coupled.  Furthermore,  for  the  case  of  constant  heat  transfer 
and  mass  injection  rate  at  the  plate,  the  temperatures  and  mass 
fraction  distributions  at  the  wall  are  not  known  a  priori  and 
must  be  determined  from  local  heat  and  mass  balances  within  the 
porous  heat  generating  plate.  This  introduces  the  coolant 
reservoir  temperature,  Tc,  and  mass  fraction,  wc ,  which  may  be 
realistically  taken  as  constant.  In  the  region  between  the 
porous  plate  and  the  coolant  reservoir,  it  is  assumed  that  the 
temperature  gradients  are  negligible  and  that  the  coolant  is 
in  one-dimensional  flow  normal  to  the  plate.  Making  an  energy 
balance  on  an  infinitesimal  strip  of  the  porous  plate  and 
neglecting  axial  heat  conduction  along  the  plate,  leads  to 

9Tl  ■  *  *  ~  <r  N  ‘  qr  -  qt  (4.8) 
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where  qt  is  the  total  imposed  heat  flux  at  the  plate  taken  to 
be  constant  and  qr  is  the  total  radiant  heat  flux  which  is 
expressed  by 
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In  writing  the  expression  for  the  radiant  heat  transfer  rate  a 
number  of  assumptions  and  approximations  have  been  made.  The 
temperatures  of  the  walls  surrounding  the  plate,  Ts,  were  taken 
to  be  constant  and  the  same  on  both  the  top  and  bottom  sides  of 
the  plate.  The  plate  itself  was  treated  as  a  gray,  diffuse 
body  having  a  much  smaller  area  than  the  walls  surrounding  it. 

A  zonal  approximation  (10)  was  employed  in  writing  the  expres¬ 
sion  for  q...  The  temperatures  of  the  top  ('T')  and  bottom  I'B') 
surfaces  of  the  plate  are  different.  The  boundary  condition 
v,v  »  0  at  y  »  0  asserts  that  the  plate  is  impermeable  to  the 
main  stream  gas.  Making  a  mass  balance  on  an  infinitesimal 
control  volume  which  includes  both  the  plate  as  well  as  the  two 
interfaces,  under  the  above  assumption,  it  can  be  shown  that 
v*  and  ww  are  related  by  the  following  equation 


(wc  -  ww)pwvw  -  - 


(4.10) 


The  equation  of  state  is  given  by  the  ideal  gas  law 
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who  rc 


p  =*  pRT 

R  =  R/M  and  M  =  I/[(w,/M,)  +  (w2/M2)] 


(4.11) 


lor  a  perfect  gas  the  specific  heat  is  independent  of  tempera¬ 
ture.  The  mixing  rules  for  specific  heat,  viscosity,  thermal 
conductivity  and  diffusion  coefficient  are  available  in  the 
literature  (44).  The  variation  with  temperature  of  viscosity, 
thermal  conductivity  and  diffusion  coefficient  was  approximated 
hy  the  power  law  relations 

v/v„  =  (T/TJa,  =  (T/TJb,  D/D^  =  (T/TJC  (4.12) 

where  a,  b  and  c  are  empirical  constants  determined  from  experi¬ 
mental  data  over  the  temperature  range  of  interest. 

4.1.3  Justification  for  Decoupling  Convection  from 
Radiation 

Problems  similar  to  the  one  posed  in  the  preceding  section 
have  been  considered  by  numerous  investigators  and  extensive 
results  have  been  reported  for  both  gray  and  nongray  gases  as 
discussed  previously.  Similarity  solutions  are  not  possible 
and  the  presence  of  radiation  has  introduced  the  nonlinear 
integral  terms  in  the  energy  equation.  Approximate,  iterative 
or  finite  difference  approaches  are  some  possible  methods  for 
solving  the  equations. 

The  results  of  nongray  gas  calculations  have  shown  (45] 
that  the  use  of  gray  absorption  coefficient  in  the  analysis 
overestimates  the  interaction  of  convection  and  radiation  on 
temperature  distribution  and  wall  heat  flux.  Hence,  the  use 
of  the  gray  approximation  should  also  overestimate  the  inter¬ 
action  in  the  present  problem  in  assessing  the  coupling  between 
convection  and  radiation.  To  be  specific,  consider  C02  gas  to 
he  gray.  Analyses  have  shown  (10]  that  in  addition  to  the 
classical  dimensionless  groups  there  arise  several  additional 
parameters  which  govern  the  flow  and  temperature  fields.  The 
parameters  are  the  wall  emissivity,  and  the  dimensionless 
numbers : 


N  -  k.>c./4«Ti 
5  •  2oTiic.x/p.u.cpi> 

The  parameter  N  represents  the  relative  importance  of  conduction 
versus  radiation,  whereas  C  is  a  measure  of  the  relative  import¬ 
ance  of  the  radiation  flux  to  the  enthalpy  flux.  The  optical 
thickness  based  on  the  boundary  layer  thickness  is  estimated  to 
be  ■  *6  ■  C(PrNC)*'2  where  the  constant  C  is  approximately 
in  the  range  of  5  <  C  <  10  depending  on  PwVw/P«u»-  For  the 
conditions  expected  in  the  experimental  phase  of  the  program  of 
interest  (M  <  0.3,  Re  <  S  *  10s,  Pr  <  1 ,  Tp  <  1500°K  and 
( pwvw)/ (P^u.)  <  0.05)  it  can  readily  be  shown,  using  the  data 


SO 


for  N  [10],  that  the  boundary  layer  can  be  definitely  considered 
as  thin,  xg  <<  1.  The  results  of  extensive  calculations  by 
numerous  investigators  [14,  15,  45,  46]  for  an  optically  thin 
boundary  layer  with  negligible  viscous  dissipation,  hot  wall 
and  no  radiation  incident  on  the  boundary  layer  from  some 
external  source  indicate  that  the  effect  of  radiation  on  tem¬ 
perature  distribution  is  negligible.  Hence,  for  the  present 
conditions  it  is  justifiable  to  compute  the  temperature  distri¬ 
bution  by  decoupling  the  radiation  from  the  other  modes  of 
energy  transfer  and  neglecting  the  term  -8F/9y  in  the  energy 
Eq.  (4.3). 

Additional  justification  for  neglecting  the  radiation  term 
in  the  energy  equation  is  given  in  Figs.  15  and  16  which  show 
that  for  near-stagnant  C02  gas  the  temperature  profiles  differ 
only  by  about  II  from  those  when  radiation  was  neglected. 
Finally,  the  results  of  Schimmel  et  al.  [26]  who  have  studied 
heat  transfer  in  a  plane  layer  of  stagnant  C02  gas  bounded  by 
two  opaque  walls  show  that  the  interaction  is  very  small  for 
wall  emissivities  greater  than  about  0.2,  even  when  energy 
transfer  by  convection  is  absent. 

4.1.4  Integral  Forms  of  the  Conservation  Equations 


The  mathematical  difficulties  of  an  exact  solution  can  he 
circumvented  by  an  approximate  analysis  which  simplifies  the 
mathematical  manipulations.  The  justification  for  doing  this 
is  that  in  cases  where  exact  solutions  are  available,  they 
agree  with  satisfactory  accuracy  with  the  solutions  obtained 
by  the  approximate  method. 


Multiplying  the  continuity  Eq.  (4.1)  by  u,  then  adding  to 
the  momentum  Eq.  (4.2),  integrating  the  latter  with  respect  to 
y  between  the  limits  y  •  0  and  y  ■  6V»  and  by  utilizing  the 
boundary  conditions,  Eq.  (4.S)  leads  to 
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Similarly,  integrating  the  energy  Eq.  (4.3)  over  the  thermal 
boundary  layer  thickness  between  y  *  0  and  y  *  6t,  satisfying 
the  boundary  conditions  Eq.  (4.6)  and  noting  that  due  to  the 
dependence  of  cp  on  composition,  the  enthalpy  differential  is 
given  by  v 

dh  *  c  dT  ♦  c  (1  -  c  )Tdw  (4.14) 
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with  *  Cp  /Cp  ,  the  integral  relation  for  the  energy  equa¬ 

tion  can  be  written  in  the  form 
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I  inally,  integrating  the  species  hq.  (4.4)  over  the  dif- 
fusion  boundary  layer  thickness  between  y  =  0  and  y  =  6m,  and 
by  satisfying  the  boundary  conditions,  l;q .  (4.7),  the  integral 
form  of  the  species  equation  is  obtained 
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(4.10) 


Since  the  density  varies  appreciably  over  the  boundary 
layer  thickness,  Kqs.  (4.13),  (4.14)  and  (4.15)  can  be  handled 
more  conveniently  by  introducing  the  Dorodnitsyn  transformation 
[47).  A  new  variable  t  is  defined  such  that  at  a  given  value 


of  x 


dt  =  (p/pjdy 


(4.17) 


The  velocity  boundary  layer  thickness,  for  example,  is  defined 
by  the  inverse  transformation 

f6 

5v  =  (P„/P)dt  (4.18) 

Jo 

If  Lq.  (4.17)  and  (4.18)  are  applied  Eqs.  (4.13)  to  (4.15)  can 
be  written  in  dimensionless  form  as 


A  37*  (AIV  "  *A  +  (PwPw/P.Pj  w 


(4.19) 


37*  (AFj)  *  *A(1  *  Tw/T-)  *  Pr^(Dwkw/p-k»)  [^(T/T«)]W 
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A  37*  (AF>)  "  *A(w-  '  "w}  *  5T*(0wDw/p-d-)[Ik]w  (4-21) 
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We  note  that  the  velocity  boundary  layer  thickness  6V  is 
given  by  the  inverse  transformation,  see  Eq.  (4.17), 

6y  -  [^(o./ojdt  -  6  f  (M.T/ffTJdn-  «f  ((1-w)  ♦  wfMj/MjJHT/T.ldn 
4  *  o  4 

(4.23) 
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Hence,  after  the  temperature  distribution  has  been  determined 
the  relationship  between  the  velocity  boundary  layer  thickness 
<5V  and  the  transformed  boundary  layer  thickness  6  can  be  found. 
Similar  relations  can  be  derived  between  and  5  as  well  as 
6m  and  6.  For  the  flow  treated  here,  however,  the  three  thick¬ 
nesses  (6V,  6t,  and  6m)  cannot  be  expected  to  differ  greatly. 
Moreover,  in  the  integral  method  there  appear  to  be  good  mathe¬ 
matical  reasons  [48j  for  adopting  the  point  of  view  of  one 
thickness  within  which  all  fluid  properties  change  from  their 
wall  values  to  their  free  stream  values. 

4.1.5  Solution  of  Conservation  Integral  Equations 

As  is  customary  in  the  solution  of  Karman- Pohlhausen  inte¬ 
gral  boundary  layer  equations  we  approximate  the  velocity,  tem¬ 
perature  and  mass  fraction  profiles  by  nth  order  polynomials: 
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The  coefficients  a j ,  b^,  and  cj  are  determined  so  that  the 
boundary  conditions  in  the  k*-n  plane  corresponding  to  those 
[given  by  Eqs.  (4. 5-4. 7)]  in  the  x-y  plane  are  satisfied  for 
any  value  of  6  and  for  any  value  of  Tw  and  ww.  Previous 
analyses  have  shown  that  sufficient  accuracy  for  the  purpose 
of  the  present  investigation  can  be  realized  if  fourth  order 
polynomial?  arc  assumed. 


The  boundary  conditions  on  u,  T  and  w  in  the  x*-n  plane 
corresponding  to  those  in  the  x-y  plane  given  by  Eqs.  (4.5), 
(4.6)  and  (4.7)  are 
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The  second  boundary  condition  at  the  wall  (n  =  0)  follows  from 
the  constraint  that  the  shear  stress  tw  ,  the  heat  flux  qw  and 
the  mass  flux  i%  are  constant  in  the  vicinity  of  the  wall. 

With  Eqs.  (4 . 27) - (4 . 29)  the  coefficients  for  the  velocity, 
temperature  and  mass  fraction  profiles  can  be  evaluated  and 
they  are 
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The  temperature  and  mass  fraction  gradients  at  the  wall  are 
found  from  Eqs.  (4.9)  and  (4.10).  In  terms  of  the  variables 
here  there  result 
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(4.35) 


With  Eqs.  (4.30),  (4.31)  and  (4.32)  the  F  integrals  can  be 
evaluated,  and  we  obtain 


(4.38) 

Since  the  mass  injection  rate  .  ong  the  plate  is  uniform, 

Pwvw  =  constant,  both  T,.  and  ww  will  vary  in  the  x  direction. 
Substituting  Eqs.  f 4 . 3 0 J  ,  (4.31)  and  (4.35)  into  Eqs.  (4.19), 
(4.20)  and  (4.21),  leads  respectively  to 

a?r  O', A)  -  ♦*  *  (.wUw/p.Uj(2Xv/(l*Xv))  (4.39) 
jpr  ( F2A)  ■  4U  -  (TW/TJ  *  (c  /cp)|(Tw/TJ  -  r|}  (4.40) 

t1^F  (F,A)  -  ♦(».  -  wc)  (4.41) 

which  define  A,  Tw  and  ww  as  functions  of  x.  For  an  arbitrary 
variation  of  physical  properties  with  temperature  and  the  fact 
that  the  parameter  r  is  a  nonlinear  function  of  T\.  no  closed 
form  solutions  of  the  differential  equations  (4.59)  through 
(4.41)  are  possible.  Numerical  integration  of  the  equations  is 
quite  tedious  and  the  effects  of  the  various  parameters  cannot 
be  readily  determined.  For  this  reason  some  justifiable 
assumptions  and  approximations  will  be  made  ir.  order  to  eliminate 
the  simultaneous  solution  of  the  equations.  There  are  3lso  seme 
difficulties  associated  with  the  boundary  conditions  which  will 
be  discussed  before  the  solutions  can  be  obtained. 

The  radiant  heat  flux  qr  appearing  in  the  parameter  r 
depends  on  x  since  Tw  varies  along  the  plate;  however,  for  the 
conditions  of  interest  qr  is  considerably  smaller  than  q*.  For 
this  reason  and  the  fact  that  Eq.  (4.9)  is  only  approximate  in 
the  first  place  q^  will  he  taken  as  constant.  The  assumption 
of  constant  qr  and  hence  of  T  is  justified  on  the  basis  that 
the  gray,  diffuse,  lumped  parameter  approximations  made  in 
deriving  liq.  (4.9)  and  the  uncertainty  in  the  values  of  the 
emissivity  arc  more  serious  than  the  assunption  of  qr  =  const 


at  the  plate.  The  approximation  can  be  improved  by  basing  q 
on  some  mean  plate  temperature.  The  physical  property  variation 
with  temperature  is  approximated  by 

v/vm  -  k/k.  =  T/T.,  D/D^  =  (T/TJS,  cpj/cp  =  1  (4.42) 

l-'or  the  temperature  range  of  interest  the  values  of  exponents 
a  =  b  =  1  and  c  =  2  are  slightly  high.  A  somewhat  better 
approximation  would  be  a  =  b  =  3/4  and  D  =  3/2.  Unfortunately, 
for  this  latter  case  an  analytical  solution  would  not  be  possi¬ 
ble.  It  is  further  assumed  that  in  the  expressions  for  Xv,  X^ 
and  Xm,  Eq.  (4.33),  the  following  approximations  can  be  made: 
r(Uoo/uw)  -  T,  r(k„/kw)  =  r  and  rfDoo/Dw)  ~  r.  Again,  in  view  of 
the  assumptions  already  made  this  appears  to  be  a  reasonable 
compromise . 

At  the  leading  edge  x  =  0  and  at  x  =  x0  a  discontinuity  in 
the  solution  exists  because  the  boundary  conditions  are  not 
continuous.  Indeed,  it  is  possible  that  in  the  immediate 
neighborhood  thereof  the  boundary  layer  equations  themselves 
may  not  be  valid  since  the  boundary  layer  hypothesis  nrecludes 
such  abrupt  change  in  the  flow  variables  in  the  x-d\rection. 

In  the  analysis  the  discontinuity  has  been  handled  by  assuming 
that  the  boundary  layer  equations  are  valid.  This  approach  can 
be  expected  to  be  valid  provided  the  discontinuities  of  the 
local  velocity,  temperature  and  mass  fraction  within  the  bound¬ 
ary  layer  are  relatively  small  [47] . 


The  solutions  of  the  equations  corresponding  to  the  start¬ 
ing  conditions  illustrated  in  Fig.  24a,  for  example,  will  be 
presented.  Substituting  Eqs.  (4.32)  ar.d  (4.38)  into  Eq.  (4.39) 
and  making  use  of  the  above  assumptions  results  in 
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Integration  of  Eqs.  (4.40)  and  (4.41)  yields  respectively 
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In  Eqs.  (4.43)  through  (4.45)  50(=  x*<|>2)  is  the  dimensionless 
distance  along  the  plate  at  which  the  injection  of  the  coolant 
begins,  and  Xv  0  is  the  corresponding  dimensionless  thickness 
of  the  boundary  layer.  The  parameter  C0  can  be  determined  by 
going  to  the  limit  4>  -*•  0  [49].  More  simply,  the  limiting  case 
of  <j>  -*■  0  corresponds  identically  to  that  of  an  impermeable  wall 
for  which  the  boundary  layer  thickness  is  [50] 

6(x0)/x0  =  5.67/^AT  (4.46) 

Therefore  the  beginning  of  the  injection  region  is  given  by  the 
following  expression 

c.  =  /VJ  -  (4.47) 

Equations  (4.43),  (4.44)  and  (4.45)  together  with  the  equations 
for  the  F's,  Eq.  (4.36)  through  Eq.  (4.40),  complete  the  solu¬ 
tion  for  i|»,  Tw/T^,  and  ww  as  functions  of  x,  i.e.  ,  =  £(i]0, 

Tw/T<n  *  Tw  (C  ,<J0 /To,,  and  ww  =  w„(?,  T^/To,) .  The  velocity, 

temperature  and  mass  fraction  distributions  can  be  computed 
from  Eqs.  (4.30)  -  (4.32).  The  solutions  corresponding  to  the 
starting  conditions  illustrated  in  Fig.  24b  can  be  given  in  a 
similar  manner. 

4.1.6  Prediction  of  Attenuation 

As  already  discussed,  difficulties  were  encountered  in 
making  a  direct  comparison  between  the  predicted  and  measured 
attenuation  when  using  the  total  band  absorptance  model  since 
the  band  width  parameters  could  not  be  adjusted  to  correspond 
with  those  of  the  wide  band  path  filter.  For  this  reason  the 
spectral  absorption  coefficient  instead  of  the  band  absorptance 
is  used  to  predict  the  band  attenuation  as  discussed  below. 

Since  the  intensity  of  the  external  radiation,  Ie\,  is  much 
smaller  than  the  intensity  of  radiation  emitted  from  the  plate, 
EXIbx(’rp)*  the  second  term  on  the  right-hand-side  of  Eq.  (2.2) 
can  be  neglected.  With  this  approximation  the  radiation  leaving 
a  diffusely  emitting  and  reflecting  plate  reduces  to 

r 6  -/yK,dy/cos6 

Iu(°.8)  -  exIbx(Tp)  *  cx  j  tbx(y)e  Kxdy/cose  (4.48) 

0 

and  for  a  specularly  reflecting  plate  having  both  directional 
dependent  emissivity  and  reflectivity  the  spectral  intensity 
leaving  becomes 
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With  IU(O,0)  determined,  the  spectral  normalized  intensity 
lA*(<$,0)  and  the  band  intensity  I*(6,0)  can  be  predicted  from 
Eqs.  (2.4)  and  (2.5),  respectively,  if  the  spectral  absorption 
coefficient  is  known. 

Since  the  integration  is  along  a  nonisothermal  and  non- 
homogeneous  path,  the  spectral  absorption  coefficient  k\  must, 
be  known  as  a  function  of  temperature  and  composition.  The 
values  of  were  therefore  taken  from  the  computations  of 
Malkmus  [51]  who  presented  spectral  emissivity  data  as  a  function 
of  temperature  for  the  weak  and  strong  line  limits.  Briefly,  he 
assumed  a  harmonic  oscillator  for  the  C02  molecule  and  a  random 
Elsasser  array  for  the  line  structure.  The  assumptions  and  the 
models  are  more  completely  discussed  elsewhere  [2,  52].  In 
general,  the  weak  and  strong  line  limits  are  valid  when  the 
emissivity  (or  absorptivity)  of  the  band  center  is  well  below 
and  near  unity,  respectively.  The  use  of  the  results  of  Malkmus 
is  somewhat  arbitrary.  However,  due  to  the  lack  of  experimental 
data  in  the  temperature  range  of  interest,  the  fact  that  other 
computations  [53-54]  differ  little  from  his  results  and  they 
compare  well  with  the  available  experimental  data  [55]  the  values 
of  <x  should  be  quite  adequate. 

4. 2  Experiment 

4.2.1  Description  of  Experimental  Facility 

The  principal  equipment  used  in  the  experimental  work  of 
this  research  program  was  a  continuously  operating,  open  circuit, 
subsonic  wind  tunnel  with  a  design  capacity  of  approximately 
M  «  0.2  at  a  volumetric  flow  rate  of  12,000  cubic  feet  per 
minute.  A  sketch  of  the  subsonic  wind  tunnel  is  shown  in  Fig. 

25.  The  wind  tunnel  consists  of  a  blower,  numerous  flow 
straighteners ,  a  test  section,  an  intake  bypass  valve  (butter¬ 
fly  gate)  and  a  blower  control  gate. 

The  blower  comprises  an  air  compressor  and  an  electric 
motor.  The  compressor  is  a  Spencer  Turbo-Compressor.  It  is  a 
high  capacity,  low  pressure  differential  fan.  The  electric 
motor  is  of  the  line  start  induction  type  and  is  manufactured 
by  Allis  Chalmers.  This  motor  is  designed  for  300  HP  at  1760 
RPM  while  drawing  70  amperes  of  3  phase  60  cycle  current  at 
2300  volts. 

The  flow  straighteners  are  made  of  cardboard  tubes  that 
have  been  bound  together  and  placed  at  the  tunnel  entrance 
section.  They  are  necessary  to  reduce  entrance  turbulence 
which  could  propagate  and  produce  undesirable  effects  in  the 
test  section. 

The  wind  tunnel  test  section  is  located  on  the  intake  side 
of  the  tunnel  and  is  one  foot  square  (the  same  dimensions  as 
the  intake  duct)  and  has  three  sides  constructed  from  plexi¬ 
glass  to  facilitate  visual  and  photographic  observation.  One 
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Figure  25.  Schematic  Diagram  of  the  Facility 


side  wall  is  made  of  sliding  plexiglass  mounted  on  gear  racks 
and  serves  as  an  accers  "o  the  test  section  as  well  as  pitot- 
static  tube  scanner  support  in  both  the  horizontal  and  vertical 
directions.  The  plexiglass  plate  in  the  top  panel  of  the  test 
section  has  a  viewing  port  for  the  radiometer  and  is  also  pro¬ 
vided  with  access  ports  for  velocity,  temperature  and  gas 
sampling  probes. 

The  intake  bypass  valve  serves  for  the  control  of  the  flow 
rate  and  consequently  the  velocity  in  the  test  section  while 
operating  the  blower  at  rated  capacity.  The  blower  control 
gate  is  used  primarily  for  bringing  the  flow  rate  up  to  a 
desired  capacity  during  the  starting  cycle.  Starting  of  the 
motor  must  be  done  at  minimum  load;  hence  the  blower  control 
gate  has  to  be  kept  closed  until  the  motor  starter  has  deacti¬ 
vated.  The  gate  is  then  slowly  opened  to  achieve  the  desired 
flow  capacity.  Figure  26  is  an  overall  photograph  of  the  test 
facility. 

4.2.2  Test  Assembly 

The  test  assembly  consists  of  the  traversing  mechanism, 
support  stand,  and  test  section.  A  photograph  of  the  test 
apparatus  is  shown  in  Fig.  27. 

The  traversing  mechanism  is  the  same  yoke-shaped  one  used 
in  the  near- stagnant  experiments  and  straddles  the  test  section. 
Both  the  radiation  detector  and  the  gas  sampling  probe  were 
mounted  on  the  yoke.  The  design  and  the  functions  of  the 
traversing  have  already  been  described  in  Section  3.2.1.  The 
traversing  mechanism  was  again  mounted  on  a  suitable  support 
stand,  see  Fig.  27. 

Two  plate  configurations  were  considered  in  the  experimental 
investigation.  The  first  plate  arrangement  was  made  of  a  one- 
inch-long  starting  length  followed  by  a  three- inch- long  porous 
heated  plate.  A  schematic  sketch  of  this  system  is  shown  in 
Fig.  24a.  The  second  configuration  comprised  a  three- inch- long 
porous  starting  length  joined  to  a  three- inch- long  solid  heated 
plate.  The  sketch  of  this  arrangement  is  shown  in  Fig.  24b. 

The  plates  were  electrically  heated  with  two  DC  arc  welders 
connected  in  parallel.  The  radiation-absorbing  gas  was  supplied 
from  a  cylinder  equipped  with  a  suitable  metering  system  which 
has  already  been  described  in  Section  3.2.3. 

4.2.3  Instrumentation  and  Calibration 

The  techniques  for  the  measurement  of  the  plate  surface 
temperature,  gas  temperature,  injection  flow  rate,  power  and 
the  radiative  flux  have  already  been  described  in  Sections 
3.2.2  and  3.2.3  and  will  not  be  repeated  here. 
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The  aspirating  probe  for  the  gas  temperature  measurements 
was  adapted  to  a  gas  chromatograph  to  determine  the  concentra¬ 
tion  of  the  injected  species.  For  this  reason,  this  probe  is 
also  referred  to  as  the  gas  sampling  probe.  It  was  mounted  on 
the  traversing  mechanism  and  could  be  accurately  positioned  in 
the  boundary  layer  with  an  adjustable  rod. 

Standard  pitot-static  probes  were  used  in  conjunction  with 
water  filled  manometers  to  determine  velocity  profiles  in  the 
test  section.  The  performance  characteristics  of  the  wind 
tunnel  were  determined  and  it  was  found  that  the  velocity  pro¬ 
files  in  the  test  section  were  sufficiently  uniform  under  all 
conditions  investigated. 

Schlieren  phonographs  were  taken  in  the  vicinity  of  the 
plates  under  various  injection  rates  and  wind-tunnel  operating 
conditions.  The  Schlieren  apparatus  provides  a  qualitative 
picture  representative  of  the  variation  of  density  throughout 
the  flow  field  of  interest  [56].  The  arrangement  of  the  partic¬ 
ular  Schlieren  apparatus  used  in  this  investigation  is  schematic¬ 
ally  illustrated  in  Fig.  25. 

The  concentration  of  carbon  dioxide  was  measured  with  a 
thermal  conductivity  vapor  phase  fractometer,  commonly  referred 
to  as  a  gas  chromatograph  [57].  It  is  similar  in  its  design 
features  to  a  Perkin-Elmer  Model  154-B  Vapor  Fractometer.  A 
schematic  diagram  of  the  gas  chromatogram  is  depicted  in  Fig. 

28.  The  concentration  measurements  were  made  by  withdrawing 
samples  of  the  gas  mixture  from  different  locations  above  the 
test  plate.  The  sampling  technique  is  illustrated  in  Fig.  29. 

The  instrument  was  calibrated  with  three  different  air-carbon 
dioxide  mixtures  of  known  composition  supplied  by  the  Matheson 
Gas  Co.  A  conventional  procedure  [57]  was  employed  to  calibrate 
the  instrument. 

4.2.4  Test  Procedure 

The  velocity,  temperature  and  concentration  fields  as  well 
as  attenuation  of  radiation  leaving  the  test  plate  were  measured 
in  a  series  of  experiments.  Depending  on  the  particular  quantity 
measured,  the  pertinent  equipment  was  properly  arranged,  aligned 
and  the  wind  tunnel  started.  As  soon  as  the  desired  test  condi¬ 
tions  were  achieved,  the  specific  measurements  were  initiated. 

For  instance,  the  structure  of  the  boundary  layer  for  vari¬ 
ous  free  stream,  injection  and  heating  conditions  was  studied 
with  a  Schlieren  apparatus.  In  case  of  adiabatic  experiments, 

C02  was  injected  into  the  free  stream  at  a  preselected  rate. 

The  development  of  concentration  boundary  layer  was  observed 
visually  on  a  frosted  glass  screen  in  the  image  plane  of  the 
camera  lens.  Upon  final  adjustments  of  the  position  of  the 
knife  edge  for  maximum  contrast  of  the  image,  Polaroid  film 
was  exposed,  developed  and  fixed  in  standard  manner. 


63 


64 


Figure  28.  Schematic  Diagram  of  the  Gas  Chromatograph 
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The  carbon  dioxide  concentration  measurements  in  the 
boundary  layer  were  made  subsequently.  The  gas  sampling  p^cbe 
which  is  similar  in  construction  to  the  Pitot  probe  with  1/16 
inch  diameter  opening  was  attached  to  the  top  panel  of  the  test 
section  and  could  >e  positioned  vertically.  An  air  diffuser 
which  is  an  integral  part  of  the  gas  chromatograph  was  activated 
causing  a  small  sample  of  gas  to  flow  from  the  boundary  layer 
through  the  instrument.  The  flow  rate  of  the  sample  gas  was 
measured  v\ich  a  rotometer,  see  Fig.  29,  and  adjusted  to  the 
desired  level  by  controlling  ':he  discharge  rate  of  the  compressed 
air  flowing  through  the  diffuser.  The  sample  was  isolated  and 
passed  with  the  carrier  gas  (He)  through  the  column  of  the  gas 
chromatograph.  Finally,  the  fractogram  was  recorded  with  a 
strip  chart  recorder  and  analyzed. 

The  gas  sampling  probe  was  also  used  as  an  aspirating  probe 
for  temperature  measurement  in  the  boundary  layer.  After  a 
desired  plate  surface  temperature  was  achieved  by  proper  adjust¬ 
ment  of  the  electrical  power  supply  (Westinghouse  DC  arc  welders) 
the  temperature  and  velocity  profiles  in  the  boundary  layer  were 
measured  at  a  given  distance  from  the  leading  edge.  Unfortunately, 
the  plate  surface  temperature  was  limited  by  the  inadequate  power 
available  at  the  wind  tunnel  which  prevented  the  study  of  attenu¬ 
ation  over  a  wide  enough  range  at  experimental  conditions. 

Upon  completion  of  the  desired  velocity  and  temperature 
measurements  the  probe  was  removed  and  the  radiation  detector 
described  earlier  was  positioned  above  the  plate  without  inter¬ 
fering  with  the  flow  field.  As  a  consequence  of  the  construc¬ 
tion  features  of  the  wind  tunnel  only  near  normal  (6  <  20°) 
radiation  flux  measurements  were  feasible.  Once  it  was  decided 
upon  plate  configuration,  free  stream,  injection  and  surface 
temoerature  conditions,  attenuation  measurements  were  first  made 
with  a  nonabsorbing  gas  (air)  and  then  repeated  under  otherwise 
identical  conditions  by  injecting  an  absorbing  gas  (C02). 

4 . 3  Results  and  Discussion 

4.3.1  Analytical  Results 

Since  there  is  an  even  large;  number  of  independent  param¬ 
eters  than  for  the  near-stagnant  gas  layer  [i.e.,  Rex,  ♦>  C,  Pr, 

Sc,  r,  c0,  Tw,  Tc,  Ts,  wc,  plate  arrangement,  radiation  charac¬ 
teristics  of  the  plate,  incident  radiation,  radiation  character¬ 
istics  of  the  gas,  physical  properties  of  the  fluids  (constants 
a,  b,  and  c) ,  etc.],  it  is  impractical  to  cover  a  complete  range 
of  parameters  which  may  be  of  possible  interest.  For  this 
reason  only  sample  results  are  presented  for  the  diffusely  emit¬ 
ting  and  reflecting  plate  in  the  absence  of  external  radiation 
with  air  as  the  free  stream  fluid  (Pr  *  0.70)  snd  carbon  dioxide 
(Sc  =  0.5)  as  the  iniected  fluid.  Tne  effects  of  plate  emis- 
sivlty  and  emergent  direction  have  already  been  discussed  at 
some  length  in  Section  3.1.3  and  therefore  will  not  be  repeated. 
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In  the  ensuing  discussion  only  the  effects  of  flow  and  injection 
variables  upon  the  attenuation  for  given  heating  conditions  will 
be  emphasized  since  the  mass  fraction  and  temperature  of  the 
absorbing  gas  in  the  boundary  layer  are  functions  of  these 
parameters . 

The  solutions  as  calculated  from  Eqs .  (4.44),  (4.45), 

(4.46)  and  (4.47)  are  shown  in  Figs.  30  and  31  for  starting 
length  parameters  of  Co  =  0  and  £0  =  5,  respectively.  In  Figs. 
30a  and  31a  the  dimensionless  parameter  if)  which  is  a  blowing 
Reynolds  number  based  on  6  as  the  characteristic  dimension  is 
plotted  against  the  parameter  z,  (=  Rex<t>2).  Inspection  of  the 
figures  reveals  that  the  relation  between  \p  and  C  is  linear 
except  for  C  <  1.  The  variation  of  the  parameters  Av,  A^  and 
Am  with  C  is  illustrated  in  Figs.  30b,  30c  and  30d,  respectively. 
These  parameters  characterize  the  velocity,  temperature  and 
concentration  boundary  layer  thicknesses,  see  Eq.  (4.34).  It 
is  seen  from  the  figures  that  the  dimensionless  boundary  layer 
thicknesses  Av,  At  and  Ar  all  increase  with  the  increase  in  the 
dimensionless  heat  flux  parameter  T.  We  note  that  T  =  1  corre¬ 
sponds  to  adiabatic  flow  over  a  plate  which  is  maintained  at 
the  same  temperature  as  the  main  free  stream  gas.  When  the 
injection  of  fluid  begins  at  a  certain  distance  from  the  leading 
edge  of  a  plate  the  boundary  layer  thickness  at  the  beginning 
of  the  porous  section  is  taken  to  be  that  at  the  end  of  the 
impermeable  section.  This  explains  the  difference  between  the 
results  in  Figs.  30  and  31. 

The  variation  of  the  temperature  and  the  mass  fraction  of 
the  injected  species  along  the  wall  are  illustrated  in  Figs. 

32a  and  32b,  respectively.  Both  the  wall  temperature  and  the 
wall  mass  fraction  of  C02  increase  with  £.  The  mass  fraction 
at  the  wall,  ww,  for  T  ■  4  was  very  close  to  that  for  T  -  2  and 
separate  curves  could  not  be  drawn.  As  expected  the  effect  of 
T  is  much  greater  on  Tw  than  on  ww.  As  C  becomes  large, 

Tw/Tv  ■*  r  and  ww  1.  This  limiting  behavior  is  expected  on 
physical  grounds. 

The  variation  of  the  temperature  profiles  in  the  boundary 
layer  is  illustrated  in  Fig.  33.  The  dimensionless  velocity 
(u/u„,)  and  mass  fraction  ( (w-ww)/ (wB-ww)  ]  profiles  in  this  type 
of  plot  are  within  a  fraction  of  a  percent  from  those  of  the 
temperature  for  corresponding  values  of  $  and  Rex,  and  hence 
separate  curves  are  not  presented.  The  abcissa  yW6t  instead 
of  y/6  was  purposely  chosen  so  that  the  profiles  for  different 
values  of  parameter  $  would  be  separated.  The  results  show 
that  the  boundary  layer  thickness  increases  with  ♦  and  r,  sec 
Figs.  33a  and  33b.  Even  though  direct  comparison  is  not  possi¬ 
ble  because  of  different  boundary  conditions  and  temperature 
dependence  of  physical  properties  considered,  the  trends  in  the 
results  presented  in  the  figure  agree  with  those  of  previous 
investigators  (49,  58,  59).  It  is  noted  that  because  of  injec¬ 
tion  of  the  fluid  along  the  plate,  inflection  points  appear  in 
the  temperature  (velocity  and  mass  fraction)  profiles.  It  is 
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Figure  31. 


Variation  of  the  Parameters  Xy, 
C-;0  for  Different  r ;  »  5.0,  Pr 


Xt  and  X_  with 
-  0.7,  Sc  -  0.5 
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Figure  33.  Temperature  Distribution  in  the  Boundary  Layer  for 
Different  C0  =  0,  Pr  *  0.7,  Sc  =  0.5 

found  that  the  larger  the  value  of  0  the  further  the  inflection 
point  moves  outwards  from  the  plate.  The  instability  of  the 
laminar  boundary  layer  may  be  interpreted  from  the  inflection 
points  occurring  in  these  curves. 

The  effect  of  injection  parameter,  0,  and  the  temperature 
(parameter  r)  on  the  attenuation  of  radiation  in  the  4 . 3y  hand 
of  C02  is  illustrated  in  Figs.  34a  and  34b  based  on  the  weak- 
line  and  strong-line  approximations,  respectively.  The  spectral 
attenuation  was  computed  by  numerically  integrating  Eq.  (2.5J 
and  using  the  spectral  absorption  coefficients  from  the  weak- 
and  strong-line  limits  as  discussed  in  Section  4.1.6.  The 
attenuation  based  on  the  weak-line  approximation  is  greater  at 
all  wavelengths  than  the  strong- line  approximation  which  is 
consistent  with  the  definition  of  (SI).  The  parameters  r  *  3 
and  r  *  4  correspond  to  the  plate  surface  temperature  of  about 
Tw  ■  1500°R  and  Tw  *  2000°R,  respectively.  It  is  seen  that  in 
this  temperature  level  range  the  broadening  of  the  band  by 
temperature  is  small.  The  band  attenuation  is  larger  for  r  *  4 
than  for  r  *  3.  The  effect  of  the  injection  parameter,  <fr,  on 
the  attenuation  is  as  expected.  With  the  increase  in  mass  of 
CO?  in  the  boundary  layer  the  attenuation  also  increases. 

The  dependence  of  the  integrated  band  intensity  1*  on  the 
parameter  c  is  shown  in  Figs.  35a,  35b,  35c,  and  35d  for 
Rex  *  2.5  x  10"  with  T  *  3  and  4  as  well  as  for  Rex  a  H)  x  l|)" 
with  r  ■  3  and  4.  The  effects  of  direction  0  on  attenuation 
have  already  been  discussed  at  some  length  in  Section  3.3.1  ami 
the  results  are  given  therefore  only  for  the  normal  (0  =  0) 
direction.  This  of  course  is  the  most  critical  since  the 
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Figure  34.  Effect  of  4  and  T  on  the  Normalized  Emergent 

Spectral  Intensity  with  C02  Injection;  x0  =  0, 

Rex  *  10  x  10\  *  530°R,  e  *  1 

attenuation  is  minimum.  It  is  seen  from  the  figure  that  for 
e  ■  1  and  small  c  the  attenuation  increases  (weak- line  approxi¬ 
mation)  rapidly  with  c  but  then  it  levels  off  and  the  additional 
mass  of  CO 2  in  the  boundary  layer  becomes  rather  ineffective  in 
attenuating  radiation  as  was  the  case  for  the  near-stagnant  layer. 
The  strong- line  approximation  predicts,  on  the  other  hand,  a  more 
gradual  increase  in  I*  with  c.  It  is  noted  that  for  e  ■  0.5 
there  is  a  maximum  in  I*  with  £•  This  is  due  to  the  emission 
as  well  as  reflection  by  the  plate  of  the  radiation  emitted  by 
CO,  in  the  boundary  layer  and  is  consistent  with  the  definition 
of  I*,  see  Eqs.  (2.1),  (4.48)  and  (2.5).  Comparison  of  results 
presented  in  Figs.  35a  and  35b  with  those  given  in  Figs.  35c 
and  35d  show  that  for  a  given  c,  radiation  is  attenuated  more 
effectively  for  lower  Rex.  That  is,  for  smaller  free  stream 
velocities,  smaller  injection  is  needed  to  achieve  the  same 
attenuation.  The  attenuation  is  larger  for  T  -  4  than  for  T  *  3. 
The  effects  of  plate  emissivity  on  I*  agree  with  those  for  the 
near- stagnant  layer;  however,  the  influence  of  e  on  I*  is  stronger 
than  that  illustrated  in  Fig.  18  but  the  trends  are  still  the 
same. 
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Figure  35.  Effect  of  Injection  on  the  Normalized  Emergent 
Band  Intensity;  x0  =  0,  T®  =  530°R 

4.3.2  Experimental  Results 

The  development  of  concentration  boundary  layers  due  to 
mass  injection  through  a  cold  porous  plate  is  illustrated  in 
Figs.  36  and  37.  The  effect  of  injection  mass  flux  pwvw  with 
Uoo  ■  50  ft/sec  is  depicted  in  Figs.  36a  and  36b.  As  expected 
the  boundary  layer  thickness  decreases  as  the  free  stream 
velocity  is  increased,  see  Figs.  36a  and  36c  and  Figs.  37a  and 
37b.  The  photographs  in  Fig.  37  show  the  increase  of  the  bound 
ary  layer  thickness  with  decrease  in  free  stream  velocity  for 
a  constant  injection  mass  flux. 
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(b) 


p  v  =  .0067  lb/£t2sec 
w  w 

uro  =  50  ft/sec 


Figure  36.  Schlerien  Photographs  of  Boundary  Layer  Development 
with  C02  Injection  through  a  Cold  Porous  Plate 
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(a) 


p  v  =  .0167  lb/ft2sec 
Kw  w 

Uoo  =  100  £t/sec 


(b) 


p  v  =  .0167  lb/ft2 sec 
w  w 

Uoo =  100  ft/sec 


p  v  »  .0167  lb/ft2sec 
w  w 

u«o  ■  25  ft/sec 


Figure  37.  Schlerien  Photographs  Illustrating  the  Effect  of 

Free  Stream  Velocity  on  Boundary  Layer  Development 
with  C02  Injection  through  a  Cold  Porous  Plate 


Figure  38.  Measured  Velocity  Profile  in  the  Boundary  Layer 
with  C02  Injection  through  a  Heated  Porous  Plate; 
xQ  ■  1  in,  x  ■  2-1/2  in,  Rex  ■  5.75  *  10\  T„  « 
530°R,  Tw  *  1460#R 

Sample  results  of  the  velocity,  temperature  and  mass  frac¬ 
tion  profiles  in  the  boundary  layer  are  displayed  in  Figs.  38, 

39  and  40,  respectively.  Because  of  the  probe  dimensions  and 
slight  vibrations  of  the  wind  tunnel  to  which  the  probe  was 
attached  the  closest  that  the  velocity,  temperature  and  concen¬ 
tration  could  be  measured  was  1/16  in.  from  the  plate.  There 
is  some  scatter  in  the  data  but  the  trends  are  consistent.  The 
mass  fraction  at  the  wall  ww  could  not  be  measured  and  the 
values  were  obtained  by  extrapolation.  The  temperature  profiles 
show  expected  trends,  but  the  gradients  near  the  wall  appear  to 
be  high.  The  large  near- isothermal  region  in  the  boundary  layer 
is  characteristic  of  turbulent  rather  than  laminar  flow  even 
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Figure  39.  Measured  Temperature  Profile  in  the  Boundary  Layer 
with  CO]  Injection  through  a  Heated  Porous  Plate; 
x0  -  I  in,  x  »  2-1/2  in,  Rex  •  5.75  *  10\  T.  • 
530*R,  Tw  »  1460*R 

though  Re  is  considerably  smaller  than  the  transition  Reynolds 
number.  The  boundary  layer  was  also  found  to  be  considerably 
thicker  than  expected.  Some  possible  reasons  for  these  findings 
may  be  the  following.  The  heat  exchange  effectiveness  of  the 
porous  plate,  i.e.,  the  compactness  of  the  plate,  was  not  high 
enough  to  assure  that  the  coolant  injected  into  the  stream  was 
at  the  plate  temperature.  Heating  caused  some  expansion  of  the 
plate  which  resulted  in  a  small  buckling  away  from  the  stream, 
and  the  joint  between  the  cold  solid  leading  edge  and  the  heated 
porous  plate  was  not  smooth  and  may  have  caused  turbulence. 
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Figure  40.  Measured  Mass  Fraction  Profile  in  the  Boundary 

Layer  with  C02  Injection  through  a  Heated  Porous 
Plate;  xo  *  1  in,  x  *  2-1/2  in,  T,,,  ■  530°R, 

Tw  =  1460°R 

The  normalized  band  intensities  for  the  case  of  injeccion 
through  a  porous  heated  plate  having  a  starting  length  of  x0  * 

1  in.  are  presented  in  Table  6  and  are  plotted  in  Fig.  41  as  a 
function  of  C.  Meaningful  attenuation  results  for  larger  values 
of  injection  parameter  $  could  not  be  obtained  because  the 
temperature  of  the  plate  became  too  nonuniform  over  the  area 
viewed  by  the  radiometer.  We  note  that  the  measured  attenuation 
at  a  fixed  position  along  the  plate  is  correlated  well  by  param¬ 
eter  ?.  This  of  course  is  somewhat  surprising  since  the  optical 
thickness  of  the  boundary  layer,  Eq.  (2.1),  is  not  only  a 
function  of  the  total  mass  of  the  absorbing  species  but  also  of 
the  concentration  and  temperature  distribution  in  the  boundary 
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Table  6.  Effect  of  Injected  Mass  Flux  (lbs/ft2sec)  on  Attenua¬ 
tion  (I*)  for  Porous  Heated  Plate  with  x0  =  1  in, 
x  =  2-1/2  in,  Tw  =  1460°R,  0=0° 


P  v 

w 

0.050 

0.075 

0.100 

2.88 

O 

X 

0.828 

0.813 

0.795 

5.75 

X  10“ 

0.857 

0.821 

0.811 

8.62 

x  101* 

0.864 

0.833 

0.818 

Table  7.  Effect  of  Injected  Mass  Flux  (lbs/ft2sec)  on  Attenua¬ 
tion  (I*)  for  Solid  Heated  Plate  with  x0  *  3  in, 
x  -  S  in,  Tw  ■  1460#R,  6-0° 


pvVw 

0.033 

0.072 

0.126 

0.224 

2.88 

*  lO* 

0.880 

0.800 

0.745 

0.700 

5. 75 

*  lO* 

0.862 

0.806 

0.746 

0.701 

1.0 


Figure  41.  Effect  of  C02  Injection  on  Measured  Band  Attenuation 
with  a  Heated  Porous  Plate;  x0  -  1  in,  x  -  2-1/2  in. 
T»  -  S30°R,  Tw  ■  1460#R 
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layer.  The  results  for  the  normalized  band  intensity  for  injec¬ 
tion  through  a  cold  porous  plate  joined  to  a  solid  heated  plate 
having  a  starting  length  of  x0  =  3  in.  are  presented  in  Table  7 
and  are  plotted  in  Fig.  42.  In  this  case  the  data  are  not 
correlated  as  well  against  the  parameter  £.  The  results  of 
both  Figs.  41  and  42  show  that  for  a  given  injection  parameter 
<p  the  effect  of  increasing  local  Reynolds  number  is  to  increase 
the  attenuation,  but  the  effect  is  small.  The  reason  for  this 
is  that  for  a  given  blowing  rate  the  mass  of  C02  in  the  bound¬ 
ary  layer  increases  with  Rex. 

4.3.3  Comparison  of  Analysis  and  Experiment 

A  comparison  of  predicted  and  measured  temperature  profiles 
is  shown  in  Fig.  43.  The  measured  temperature  gradients  in  the 
vicinity  of  the  porous  plate  were  higher  than  those  predicted, 
and  the  measured  boundary  layer  thickness  was  much  thicker  than 
the  calculated  one.  It  is  clear  from  the  figure  that  the  mathe¬ 
matical  equations  did  not  properly  model  the  physical  phenomenon. 
A  cause  for  higher  temperature  gradients  in  the  vicinity  of  the 
surface  is  believed  to  be  the  effectiveness  of  the  Poroloy  plate 
which  was  not  sufficiently  large  enough  resulting  in  the  coolant 
being  injected  into  the  stream  at  a  lower  temperature  than  the 
plate.  Measured  boundary  layers  were  thicker  than  those  pre¬ 
dicted  since  the  Poroloy  plate  used  was  not  compact  enough,  and 
the  C02  was  injected  into  the  main  stream  through  a  large  number 
of  small  holes  rather  than  uniformly.  The  injected  mass  flux 
pwvw  was  based  on  the  total  area  rather  than  area  of  the  pores. 
Based  on  an  unoxidized  plate  the  latter  was  estimated  to  be  about 
a  factor  of  25  times  smaller.  The  pore  area  of  an  oxidized  plate 
was  even  smaller.  The  velocity  of  the  C02  gas  emanating  from  the 
pores  based  on  the  actual  pore  area  was  estimated  to  be  high, 
i.e.,  in  the  range  from  15  to  60  percent  of  the  free  stream 
velocity.  Hence  it  is  justifiable  to  say  that  C02  was  injected 
into  the  boundary  layer  through  a  large  number  of  discrete  holes 
and  the  flow  was  not  uniform  but  formed  a  large  number  of  small 
interacting  jets.  These  jets  entered  the  stream  at  an  angle  of 
about  45°.  As  the  emerging  jet  of  hot  C02  gas  enters  the  stream 
it  dissipates,  spreads,  mixes  with  the  gas  in  the  boundary  layer 
and  its  momentum  carries  the  injected  species  further  away  from 
the  plate.  This  thickens  the  boundary  layer  and  may  even  cause 
it  to  become  completely  turbulent.  The  injection  could  have 
possibly  made  the  boundary  layer  unstable  and  increased  its 
thickness  considerably  beyond  what  is  predicted  by  the  analysis. 

The  predicted  mass  fraction  distributions  differ  appreciably 
from  the  measured  ones.  The  analytical  profiles  for  the  condi¬ 
tions  shown  in  Fig.  44  were  very  close  and  two  separate  curves 
could  not  be  drawn.  Again,  the  measured  concentration  boundary 
layer  thickness  is  much  larger  than  the  calculated  one.  This 
discrepancy  is  believed  to  be  the  same  as  for  the  thermal 
boundary  layer. 
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Figure  42.  Effect  of  CO2  Injection  through  a  Cold  Porous 

Leading  Edge  on  Measured  Band  Attenuation  with  a 
Solid  Heated  Plate;  x0  ■  3  in,  x  -  S  in,  T«  * 
S30°R,  Tw  *  1460°R 
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Figure  43.  Comparison  between  Analysis  and  Experiment  of 
;  Temperature  Profiles  in  the  Boundary  l.ayer  with 

I  CO,  Injection  through  a  Heated  Porous  Plate; 

x0  -  1  in,  x  -  2-1/2  in,  Rex  »  5.75  x  to",  T*,  = 
[  530° R,  Tw  «  1460°R 
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Figure  44.  Comparison  between  Analysis  and  Experiment  of  Mass 
Fraction  Profiles  in  the  Boundary  Layer  with  CO* 
Injection  through  a  Heated  Porous  Plate;  x0  *  1  in, 
x  -  2-1/2  in,  T„  •  530°R,  Tw  «  1460*R 
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A  comparison  of  measured  and  predicted  attenuation  is  shown 
in  Fig.  45.  The  predicted  I*  is  based  on  a  plate  emissivity  of 
0.7.  This  is  a  few  percent  higher  than  the  results  of  Fig.  23 
indicate  at  X  =  4.3y  for  a  less  porous  plate.  This  somewhat 
higher  value  is  justified  because  of  the  larger  cavity  effect 
of  the  more  porous  Poroloy  plate.  The  predicted  attenuation 
was  also  corrected  for  the  transmission  of  the  filter,  Fig.  13. 
This  correction,  however,  increased  I*  only  by  less  than  one 
percent  from  the  value  predicted  without  the  filter.  As  in  the 
case  of  the  near  stagnant  layer  and  other  investigators  [43,  60] 
the  measured  attenuation  is  slightly  larger  than  the  predicted 
but  the  trends  are  in  agreement. 

The  prediction  of  the  absorptivity  of  a  nonisothermal  gas 
has  been  the  subject  of  much  current  research  effort  [43,  60-65], 
Detailed  calculations  using  the  Curtis-Godson  [61,  62]  or  box 
model  [63]  have  been  made  to  predict  the  spectral  absorptivity 
and  the  Curtis-Godson  approximation  has  been  used  to  predict 
nonisothermal  total  band  absorptance  [43,  64,  65].  Agreement 
between  predicted  and  observed  results  ranges  from  20%  [60]  to 
less  than  10%  [62,  63,  64].  On  this  basis  the  comparison  of 
the  analytical  with  the  experimental  results  of  this  investiga¬ 
tion  are  considered  relatively  good. 

As  Fig.  34  shows,  the  center  of  the  band  is  not  yet  opaque, 
hence  the  weak-line  is  the  more  appropriate  approximation.  This 
is  in  agreement  with  the  observed  results.  There  is  uncertainty 
in  the  emissivity  as  it  was  impossible  to  measure  the  values  of 
the  spectral  emissivity  in  situ.  As  has  been  shown  an  error  of 
10%  in  the  emissivity  has  a  considerable  effect  on  the  value  of 
attenuation  predicted.  Also,  the  discrepancy  between  the  pre¬ 
dicted  and  measured  temperature  and  C02  concentration  profiles 
would  be  reflected  in  the  optical  thickness  of  the  layer  as 
well  as  Planck's  function.  The  larger  measured  than  predicted 
values  of  attenuation  I*  are  consistent  with  the  measured  con¬ 
centration  boundary  layers  being  thicker  than  those  predicted. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  section  conclusions  and  recommendations  are  pre¬ 
sented  which  are  based  on  analyses  and  somewhat  limited  measure¬ 
ments  of  preliminary  nature  conducted  in  assessing  the  effective¬ 
ness  of  an  absorbing  gas  in  attenuating  infrared  radiation 
emitted  from  hot  surfaces. 

An  apparatus  for  the  investigation  of  attenuation  of  infra¬ 
red  radiation  emitted  from  a  hot  porous  surface  through  which 
an  absorbing-emitting  gas  was  transpired  was  designed  and  con¬ 
structed.  Instrumentation  and  measurement  techniques  were 
developed  for  the  determination  of  attenuation,  flow  and  state 
variables  in  a  near  stagnant  layer  of  gas  and  boundary  layer 
flow  over  a  surface.  In  a  parallel  study  analyses  were  per¬ 
formed  to  predict  the  attenuation. 

Specifically,  carbon  dioxide  was  transpired  through  a 
porous  (Poroloy)  plate  into  a  near  stagnant  layer  of  air  and 
the  attenuation  of  radiation  measured  in  the  4.3y  band.  The 
effects  of  the  injected  mass  flow  rate,  plate  temperature,  gas 
layer  thickness  and  angle  of  observation  were  investigated  and 
it  was  found  that  C02  attenuated  between  54  and  44  percent  of 
radiation  in  the  wavelength  band  studied.  Although  the  analysis 
did  not  completely  simulate  the  experimental  conditions  it  never¬ 
theless  predicted  the  experimentally  observed  trends  correctly. 

The  velocity,  temperature  and  concentration  distributions 
in  a  boundary  layer  over  a  flat  plate  under  specified  free 
stream  and  injection  conditions  were  measured  and  were  in  poor 
agreement  with  predictions  based  on  an  approximate  integral 
method  of  solution  of  the  conservation  equations.  The  attenua¬ 
tion  of  infrared  radiation  for  some  conditions  reached  up  to 
301.  The  analytical  predictions  were  in  qualitative  agreement 
with  experimental  results. 

Based  on  the  results  obtained  for  the  near  stagnant  gas 
and  the  gas  flow  situations  it  can  be  concluded  that  in  a  speci¬ 
fied  wavelength  range  the  radiation  can  be  significantly  atten¬ 
uated,  provided  the  participating  gas  or  gas  mixture  strongly 
absorbs  radiation  in  this  region.  At  high  free  stream  velocities 
the  transpiration  of  an  absorbing  gas  into  a  boundary  layer  may 
not  be  sufficiently  effective  to  produce  the  desired  attenuation. 
In  order  to  enhance  attenuation  under  these  conditions  particulate 
matter  and/or  strongly  absorbing  gases,  such  as  metal  vapors 
could  be  employed. 
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Although  the  feasibility  of  the  attenuation  technique  has 
been  demonstrated,  a  number  of  problem  areas  remain  and  should 
be  further  explored.  In  particular  it  is  recommended  that  the 
following  work  be  done: 

1.  Transpiration  cooling  studies  should  be  extended  to 
situations  where  the  coolant  is  injected  into  the 
boundary  layer  through  discrete  holes  in  a  porous 
structure  (such  as  Poroloy)  instead  of  uniformly  to 
determine  the  fluid  flow  and  heat  transfer  charac¬ 
teristics  of  the  boundary  layer  as  well  as  of  the 
heat  exchange  effectiveness  of  the  material. 

2.  An  improved  experimental  apparatus  and  instrumenta¬ 
tion  should  be  developed  with  which  a  much  wider 
range  of  experimental  conditions,  geometries  and 
injected  fluids  can  be  investigated. 

3.  .Since  data  on  the  spectral  absorption  characteris¬ 
tics  of  gases  or  gas  mixtures  as  a  function  of 
pressure  and  temperature  are  very  limited  these  data 
should  be  acquired.  Especially  it  is  recommended 
that  the  properties  of  mixtures  of  gases  and  particu¬ 
late  matter  be  studied. 

4.  Measurement  of  the  attenuation  of  radiation  on  a 
spectral  bases  in  nonhomogeneous  media  in  which 
large  temperature  gradients  prevail  should  be 
performed. 
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